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We investigate the field theory of strings having as a target space an arbitrary discrete one- 
dimensional manifold. The existence of the continuum limit is guaranteed if the target space is 
a Dynkin diagram of a simply laced Lie algebra or its affine extension. In this case the theory 
can be mapped onto the theory of strings embedded in the infinite discrete line 7L which is the 
target space of the SOS model. On the regular lattice this mapping is known as Coulomb gas 
picture. Introducing a quantum string field ^ x {f) depending on the position x and the length I 
of the closed string, we give a formal definition of the string field theory in terms of a functional 
integral. The classical string background is found as a solution of the saddle-point equation which 
is equivalent to the loop equation we have previously considered [El]. The continuum limit exists 
in the vicinity of the singular points of this equation. We show that for given target space there 
are many ways to achieve the continuum limit; they are related to the multicritical points of the 
ensemble of surfaces without embedding. Once the classical background is known, the amplitudes 
involving propagation of strings can be evaluated by perturbative expansion around the saddle 
point of the functional integral. For example, the partition function of the noninteracting closed 
string (toroidal world sheet) is the contribution of the gaussian fluctuations of the string field. The 
vertices in the corresponding Feynman diagram technique are constructed as the loop amplitudes 
in a random matrix model with suitably chosen potential. 
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1. Introduction 

The recent progress in the theory of non-critical strings has been prepared by the 
conjecture that the functional integral of 2d gravity can be discretized as a sum over 
planar graphs . Such a discretization allowed to apply powerful methods of calculation 
borrowed from the theory of random matrices @. Some of the corresponding matrix 
models can be solved exactly after being reduced to a problem of non-interacting fermions 
in a common potential. This fermionic representation proved to be very helpful in practical 
calculations but its connection with the original model seems quite formal. It is not clear, 
for example, whether the nonperturbative phenomena in the fermion problem can be given 
an interpretation in terms of strings. 

It is therefore desirable to have a formalism in which the string excitations appear 
explicitely. An example of such a formalism applied to the string propagating in the 
one-dimensional space IR is the collective field approach of Das and Jevicki || . 

In this paper we develop a similar formalism for string theories with discrete target 
spaces. Bosonic string with discrete target spaces are interesting mainly because of their 
interpretation as theories of two-dimensional gravity coupled to matter. The (discrete) 
degrees of freedom of the matter field are the points of the target space. Because of the 
enormous symmetry of the problem, it is much easier to investigate the critical behaviour of 
statistical systems on a fluctuating surface than on the plane. The exact results obtained 
for a fluctuating surface can be translated according to ||, || to the case of a frozen 
geometry of the world sheet. 

The string models we are investigating are a generalization of the ADE strings con- 
structed in ref. 0. They are closely related to the SOS and RSOS solvable statistical mod- 
els on a regular planar lattice @. The target spaces of these models are one-dimensional 
discrete manifolds representing Dynkin diagrams of simply laced Lie algebras or their affine 
extensions. All such target spaces can be mapped onto a master target space which is the 
infinite discretized line Z. 

Since in a one-dimensional embedding space there is no room for transverse excitations, 
the states of the closed string is completely determined by the two global modes - the center 
of masses x and its length L. The string field ^(L) is an operator creating on the world 
sheet a boundary of length L and position x in the target space. The path integral for 
the correlation functions of the string field involve embeddings of the world sheet such 
that each connected component of its boundary is mapped into a single point of the target 
space. 

The string models related to the SOS model are remarkable with the possibility to 
reduce the string path integral to a simpler one involving only the global modes x and 
L, even before going to the continuum limit. As a result one obtains an effective two- 
dimensional QFT for the field fy x (L). 

The reduction of the string path integral stems from the loop, or polygon represen- 
tation of the SOS-type statistical systems @ which has been generalized to the case of 
irregular lattices in [lj. A map of the world sheet of the string in the target space can be 
represented by a collection of nonintersecting loops (domain walls) separating the domains 
of constant coordinate x. We assume that the discontinuity across a domain wall is ±1. 
Each loop configuration is characterized by its topology, the lengths of the loops, and the 
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coordinates of the connected domains. It can be imagined as a "stroboscopic picture" of a 
the evolution of one or several closed strings, the time direction on the world sheet being 
at all points orthogonal to the domain walls. The evolution of the string is decomposed 
into elementary processes represented by the domain walls and connected domains. The 
domain walls describe elementary steps in the X-space without changing the world-sheet 
geometry. On the other hand, the connected domains describe topology-changing processes 
without propagation in the embedding space. This nice factorization renders the model 
exactly solvable. 

A loop configuration is represented by a graph with lines corresponding to domain 
walls and vertices corresponding to the connected domains on the world sheet. There 
are vertices of all orders including tadpoles describing processes of creation or decay of a 
closed-string state. Each vertex is characterized also by the number H of handles of the 
connected domain on the world sheet. Vertices with H > describe processes involving 
creation and annihilation of H virtual closed strings at the same point x of the embedding 
space. 

The string path integral can be written as a sum over all loop configurations, and 
the entropy of a loop configuration factorizes to a product of weights associated with 
the domains. The weight of each domain is proportional to the entropy of a nonembedded 
random surface with the corresponding topology. The latter can be calculated, for example, 
by using the matrix model formalism. In this way the the graphs representing the loop 
configurations can be interpreted as Feynman diagrams for an ordinary field theory in the 
space of global modes (x, L) and the string path integral is reduced to a functional integral 
for the field V X (L). 

The first (and in a sense the most difficult) problem to solve is to find the classical 
string background (L) which is the solution of the classical equation of motion in 
the effective field theory. Then the propagation of the noninteracting closed string is 
determined by the gaussian fluctuations of the string field around the classical solution 
and the vertices are given by the higher order terms. The accomplishment of this program 
will be the subject of this paper. 

We are going to present in details the results announced in three short publications 



and [jlOfl , as well as a few new results as the microscopic construction of the dilute 
ADE models and a description of their multicritical points |TT . 
The paper is organized as follows. 

In section 2 we define the path integral for a string theory with discrete target space 
and establish its equivalence with an ordinary field theory for the string field ^ X (L). 
The measure in the space of world sheet geometries will be discretized by planar graphs. 
Therefore we have first to define the statistical systems related to the SOS model (known 
also as interaction-round-a-face, or IRF models 0) on an arbitrary irregular lattice S . In 
sect. 2.1 we present a generalized version of the construction of these models worked out in 
. The present version allows two nearest neighbour sites s and s' of S to have coordinates 
(heights) x and x' which either or coincide or are nearest neighbours in the X space. In 
this way we are able to achieve both the dense and dilute phases of the IRF models. The 
important special case X = TL is discussed in sect. 2.2. where we consider a nonunitary 
version of the SOS model with complex Boltzmann weights. This model renormalizes at 
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large distances onto a gaussian field with a linear term coupled to the local curvature of the 
discretized world sheet. The mapping of the IRF models onto this SOS model known as 
Coulomb gas picture, can be constructed in the same way as in the case of a lattice without 



curvature , [[L2| , [13[ . The Coulomb gas picture is based on the loop expansion which is 
explained in sect. 2.3. The partition function of the statistical system with target space 
X equals the sum over the embeddings of a world sheet with fixed geometry. In section 
2.4 we define the string path integral as a sum over the world sheet geometries and review 
its description by means of a two-component gaussian field fl4]] || || . The normalization of 



the electric and magnetic charges is fixed by the correspondence with the SOS model. In 
section 2.5 we establish the equivalence between the string path integral and the Feynman 
diagram series for the field ^ X (L) describing the dynamics of the global modes of the 
string. By means of a shift ^ X (L) — > ^ X (L) + \I/^(L) we define an improved diagram 
technique containing no tadpoles. The classical string background ^ X (L) is determined 
by the saddle-point equation. In sect. 2.6 we show that this equation is equivalent to the 
loop equation considered in 

In section 3 we study the continuum limit of the string field theory. For this purpose 
we introduce a cutoff a with dimension of length and replace the dimensionless length 
L with a renormalized length I = aL. The continuum limit can be achieved if the bare 
parameters of the string are tuned to a critical point. The critical points can be found 
as the singularities of the loop equation for the classical string background. The latter 
allows an interpretation as the partition function of a gas of nonintersecting loops on a 
fluctuating surface with the topology of a disc. In the vicinity of a critical point the volume 
of the world sheet always diverges and the different critical regimes are distinguished by 
the behaviour of the loops. The latter is gouverned by a parameter (energy) coupled to 
the total length of the loops. If the energy exceeds the entropy, then the loops remain 
small, and the possible critical regimes are the m-critical points of nonimbedded surfaces 
(sect. 3.1). This is the phase of noncritical loops. In the opposite case, considered in sect. 
3.2, the energy of the loops is not sufficient to compensate their entropy, and the loops 
form a dense critical phase filling the world sheet (sect. 3.2). Finally, the dilute phase of 
the loop gas (sect. 3.3) is achieved when the energy of the loops is tuned to its critical 
value. In the dilute phase the loops are still critical but they cover only a small fraction 
of the points of the lattice. In this phase the critical behaviour is sensitive to the choice of 
the integration measure in the space of world sheet geometries. By tuning the Boltzmann 
weights of the planar graphs we can achieve an infinite sequence of multicritical points of 
the dilute phase. The continuum limit of the loop equation (sect. 3.4) is universal ; the 
different critical regimes are classified by the possible asymptotics of the solution at small 
lengths. 

Once the string background is known, it is not difficult to find explicit expressions for 
the Feynman rules in the continuum limit. The vertices in the improved diagram technique 
(i.e., dressed by tadpoles) can be calculated as the loop correlators in a special matrix 
model constructed in sect. 4.1. The physical meaning of the coupling constant of this 
matrix model is discussed in sect. 4.2. It is coupled to a local operator which has negative 
dimension if the theory is not unitary. Its dimension is related to the fractal dimension 
of the connected domains on the world sheet. In sect. 4. 3 we give explicit expressions for 
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the planar (H = 0) vertices in the continuum limit. Since the nonplanar loop amplitudes 
in the one-matrix model are not yet known, we have only fixed the general structure of 
the nonplanar vertices. The expression for a vertex with H handles and n legs contains 
3H + n — 3 unknown coefficients. The Feynman rules simplify in the momentum space 
(p, E) where the propagator of the string field diagonalizes. The spectrum of on-shells 
states discussed in sect. 4.4 is given by the light cone in the (p, E) space. The latter 
has the geometry if a half-infinite cylinder since the discreteness of the X space leads to 
periodicity in the p-direction. Finally, in sect. 4.5 we check that the partition function 
of the noninteracting string (closed surfaces with the topology of a torus) is obtained by 
integrating over the gaussian fluctuations around the saddle point. 



2. Statistical systems on planar graphs and kinematics of strings with discrete 
target space 

2.1. The IRF height models on an arbitrary irregular lattice 



It is well known |T5] that the two-dimensional rational conformal-invariant QFT with 
central charge C < 1 are classified by the simply laced Lie algebras (i.e., these of the clas- 
sical series A n , D n , Eq, E?, E 8 ). Each of these theories can be constructed microscopically 
as a lattice statistical model whose local degrees of freedom are labeled by the points of 
the Dynkin diagram of the corresponding Lie algebra fl6|l . The statistical models associ- 
ated with dynkin diagrams of A, D : E type, or shortly ADE models represent a natural 



generalization the RSOS face models considered by Andrews, Baxter and Forrester ||17|| . 
The local degrees of freedom are attached to the sites of the lattice and interact through 
"interactions-round-a-face" (IRF) around each plaquette. 

Similarly, the extended ADE Dynkin diagrams describe conformal invariant QFT with 
C = 1 and discrete spectrum of conformal dimensions. In fact, Pasquier's construction ]T6 | 



can be applied to any one-dimensional discrete manifold X , i.e., a set of points x and links 
(xx f ) with the structure of a one-dimensional simplicial complex. The requirement that X 
is an (extended) Dynkin diagram of ADE type guarantees the existence of a scaling limit. 
Another important target spase is the infinite discretized line 2. This is the target spase 
of the model known as SOS (solid-on-solid) model. 

Before presenting the definition of the IRF model with target space X we are going 
to describe the space of excitations of this model. 

The one-dimensional discrete manifold X is defined by the set of its points x and the 
adjacency matrix C xx /, x, x' G X 

C xx ' = [the number of links connecting x and x'] (2-1) 

We assume that the target space X is represented by a nonoriented graph which implies 
symmetric adjacency matrix. 

The Hilbert space of states of the X-field in a fixed "time slice" consists of all closed 
paths in X with given length. Therefore in order to identify the ground state and the 
excited states we have to solve the problem of random motion in X. The propagation 
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kernel for a random walk on X consisting of n steps is just the n-th power of the adjacency 

(2.2) 



matrix C. Introducing the eigenvectors 

E, i 
C X x'V {p) = (3 p V {p) 



we can write the kernel K xx i = (C n ) xx r as a sum of projectors on the eigenstates 

[K n ) xx , = Y,^) nV M V (p) ( 2 - 3 ) 

v 

If X is an ADE Dynkin diagram (Fig.l), then the eigenvalues of C have the form 

f3 p = 2cos(7rp), p = m/h (2.4) 

where h is the Coxeter number and the integer m is one of the Coxeter exponents of the 
Dynkin diagram . The eigenvectors VT^ define the Fourier transform from coordinates 
x to the discrete momenta p = m/h. The ground state corresponds to the maximal 
eigenvalue of the adjacency matrix C 



(3 = (3 P0 = 2cos(7rp ), po = l/h (A,D,E) 



(2.5) 
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Fig.l : Dynkin diagrams and their Coxeter exponents 



If X is an extended ADE Dynkin diagram, then the spectrum of Coxeter exponents 
includes m = and the maximal eigenvalue of C is 



13 = 2, Po = 



(A,D,E) 



(2.6) 



1 The inverse statement reads: if all the eigenvalues of the adjacency matrix are less than 2, 
then the graph X is an either an ADE dynkin diagram or a quotent A2 n /Z2 []T8|| 
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In order to simplify the notations we shall always denote the ground state wavefunction 
by 

s* = V L) ( 2 - 7 ) 

Since the entities of the adjacency matrix are nonnegative, by the Peron-Frobenius theorem 
S x > for all xeX. 

The propagation kernel ( J2.»j| ) becomes in the limit of large proper time n a single 
projection operator to the ground state S x 

(K n ) xx ,~S x ((3) n S x . (2.8) 

It is convenient to define the wavefunctions of the excited states as 

X X { P )=V^/S X (2.9) 

They satisfy a closed algebra 

X( P )X(p>) = ^2C PP 'p"Xp" (2-10) 
p" 

and orthogonality conditions of the form 

E 5 *4>)4>') = V, T,*tp)4>) = t'ASmS*)- 1 (2.ii) 

x p 

Let us give two simplest examples. 

i) The Dynkin diagram of the algebra Ah-i is the chain of h—1 points x = 1, 2, h—1. 
This is the target space of the RSOS models The eigenvectors of C are 

Vg) = y| sin(Trpx), p = l//i, 2//i, (/i - l)//i; S * = \j\ sin(7rx//i) (2.12) 

In this case the ground state corresponds to po = 1/h. 

ii) The extended Dynkin diagram corresponding to the affine Kac-Moody algebra 
Aih-\ is the ring TLih made of 2h points. It can be represented by the set of integers 
modulo 2h. Then the connectivity matrix C xx ' reads 

^ = C'+i+C'-i ( 2 - 13 ) 

where 5^ 2h ^ is the Kronecker symbol modulo 2h. The eigenvectors of the matrix C are 

VjJ) = (2h)- 1 ' 2 exp(z7rpx),p = 0, ±i 1- (2-14) 

The target space is translational invariant and the spectrum of momenta contains the point 
p = 0. The ground state is S x = = 1/V2h. The vertex operators (order parameters) 

X( p) =e i7tpx (2.15) 



satisfy a closed algebra with fusion rules representing the momentum conservation modulo 
the period 2 of the momentum space 

<W =S^(p + p'+p") (2.16) 

The microscopic definition of the IRF models on a general irregular lattice was pre- 
sented in 0] as a direct generalization of Pasquier's construction Jl6] , the only new point 
being the explicit dependence of the Boltzmann weights on the scalar curvature. The 



effects of the curvature have been also studied in [13]. Below we present a more general 
construction including both the dense and the dilute versions of the IRF models. 

Let 5 be a two-dimensional discrete manifold such that all its faces are squares (we 
shall call them also plaquettes) . It can be constructed from an arbitrary planar graph and 
its dual one by taking the points of both graphs and adding bonds connecting all pairs of 
points which serve as extremities of mutually dual lines. A small section of such lattice is 
shown in Fig. 2 

Each configuration of the x-field on S defines a map S — > X. Sometimes the local 
values of the x-field are called heights. The map is continuous in the sense that nearest 
neighbours in X are images of nearest neighbours in S. In the original papers a stronger 
condition has been imposed, namely, that the links {ss') of S are mapped into links (xx') 
of X. Here we relax this condition by allowing two nearest neighbours in S to be mapped 
into the same point x G X. 




Fig. 2: A section of the discretized world sheet 

The statistical weight of each field configuration is a product of factors associated with 
the plaquettes and sites of the lattice S. We assume that the weights of the plaquettes 
are symmetric under cyclic permutations. This property will be of crucial importance for 
the generalization of the models to the case of an irregular lattices. The weight of the 
plaquette (g^) is defined by 
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(2.18) 



where T is the temperature of the statistical system. 

The factors associated with the sites s 6 5 reflect the response of the x-field to the 
local curvature R s which is concentrated at the sites. On an irregular lattice the local 
Gaussian curvature at the point s is defined as twice the deficit angle 



R a - tt(4 - q s ) 



(2.19) 



where q s is the coordination number at the site s (= the number of plaquettes sharing this 
site). We use the standard normalization such that 



sES 



(2.20) 



where x = 2 — 2H is the Euler characteristic of the surface S. The weight associated with 
a site s is 

O(x) = (S x ) ks/ ^ (2.21) 



The factors ( |2.21| )are trivial on a regular square lattice where all coordination numbers are 
equal to 4. 

Collecting all factors, we define the partition function of the model with target space 
X on the lattice S as the sum over all maps S —>■ X 



F ^ = E 11°^)) II n(x( Sl ),x(a 2 ),x(s3),a;(a 4 )) 



(2.22) 



/ S 2 S 3 \ 
^ s 1 s 4 



In the limit T — >• the r.h.s. of ( |2.22| ) is the partition function of the critical ADE 
models considered in |]J . In the high temperature limit T — ► oo the a;-field freezes 
to a constant and there are no long range correlations. At some critical temperature T c 
there should be a phase transition between the low-temperature (critical) and the high- 
temperature (noncritical) phase. At the point T c a new critical regime can be achieved. 
We call this phase, by reasons which will become clear later, dilute phase. The critical 
low-temperature phase will be refered also as the dense phase of the IRF models. 
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2.2. The nonrestricted SOS model with a false vacuum 
It is known 



19 , 2C 



, that the dynamics of the two-dimensional conformal theories 
can be described in terms of a Gaussian massless field with a linear term coupled to the 
local curvature. Let £ a , a = 1,2, be the local coordinates on a surface S with metric 
9ab(0'i a = 1)2 and local curvature R(£). Then the gaussian field is defined by the 
action 

A[x; 9] = \j d^^ffetl [ng g ab d a x(Od b x(0 + ia x(£)R(0] (2-23) 

where g is the coupling constant and «o is the background electric charge. 

One can map all conformal field theories with conformal anomaly C < 1 onto this 
gaussian field in the sense that the partition and correlation functions can be interpreted 
in terms of distributions of electric and magnetic charges in the gaussian theory. This 



mapping is known as Coulomb gas picture [^TJ, [^2|, [p3| , [p4[ . The electric charge e is 
carried by the vertex operator V a (^) = exp(iTiax(£)) and the charge fx at the point £ 
creates a discontinuity 2 fx along a line starting at the point £ (in fact, the discontinuity 
can be distributed among several lines forming a star with center at £). 
The (++)-component of the energy-momentum tensor is 



T(£) = —g7t 2 d^xd^x + Z7rao<9f 



(2.24) 



where we have introduced the complex variable £± = £i ± From the o.p.e. [19], [23 



2(£-£ 



f\4 



[1-6^] 
9 



+ 



A2 



nformal a 

of the vertex operators V a are 



a(2a — ot) 1 



(2.25) 



(2.26) 



±9 & ~ Z 

one reads that the conformal anomaly C of the field ( |2.23|) and the conformal dimensions 



1-6- 



a 



o , 



a (a — 2ao) 



(2.27) 



g 7 " 4g 

have the same conformal dimension and therefore can be 



The operators V a and V2 ao -a 
considered as related by charge conjugation. Their two-point function is 



(K(e)Wa(eo) = ie-c 



/ia(2a — o.)/g 



(2.28) 



where g is the coupling constant in ( |2.23|) . The charge neutrality is restored by the presence 
of the background charge — 2«o- The construction of higher correlation functions involves 
the so called screening operators that carry nonzero electric charge and have conformal 
dimension 1. By eq. (|2.27|) there are two such operators which are charge-conjugated to 
each other 

a±(a± — 2cto) = 4$; a+ + a- = 2ao (2.29) 
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The allowed charges are labeled according to the number of screening charges needed to 
neutralize the 4-point function |23] 



1 — T 1 — S 

a rs = — — a + H — a- (2.30) 

The conformal dimensions of the corresponding vertex operators form the Kac spectrum 

(ra+/2 + sa-/2) 2 - 0% 



4g 



(2.31) 



Finally, the vortex operator with magnetic charge (discontinuity) \i has conformal dimen- 
sion 

K = 9V/4-ag (2 . 32) 

The gaussian field dominance in the two-dimensional critical phenomena has its mi- 
croscopic equivalent. It happens that all target spaces with dimension not greater than 
one can be mapped onto the target space of the nonrestricted SOS model which is the 
discretized real line Z. The corresponding connectivity matrix 

C xx > = S X)X >+i +Sx,x'-i (2.33) 

has a continuous spectrum of excitations 

P = 2cos(7rp), -l<p<l (2.34) 

and its eigenvectors are plane waves 

Vgj = e i7rpx (2.35) 

The momentum space (the dual of the SOS target space) is therefore a circle with perimeter 
equal to two. 



The large distance behaviour of this SOS model is argued [|3[ |2J] to be described by 
the action ( |2.23j ) with an appropriate choice of the background charge ao and the coupling 
constant g. The background charge in the SOS model is introduced by taking an excited 
state with momentum po as a vacuum state 

S x = Vg o) = e l7rpoX (2.36) 

Then the curvature dependent-factor in the definition of the partition function can be 
written as an exponent 

n(^)^ )/47r =exp^^^)x( S )) (2.37) 

s s 
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which is the microscopic realization of the curvature-dependent term in ( 2.23 ). Knowing 
that the global curvature does not renormalize, we conclude that 



cto — Po + even integer 



(2.38) 



Since the charge in the SOS model is determined up to an even integer, we can assume 
that ao = po- 

Now it remains to fix the coupling constant g in the gaussian theory. 

First we observe that the spectrum of allowed momenta in the target space of the SOS 
model is p = mpo, m =integer. The operators ( |2.35| ) corresponding to momenta outside 
this spectrum have vanishing correlators because the electric neutrality is not fulfilled. On 
the other hand, due to the discreteness of the target space of the SOS model all charges of 
the form a ± 2n in the gaussian model are indistinguishible in the SOS model. Therefore, 
the spectrum of allowed electric charges in the gaussian model is 



a = rp + 2n; r,n G Z 



(2.39) 



The way the coupling constant g depends on the background charge ao is fixed by the 
compatibility of (|2~2^|) and ( |2~39D 



a = g-l; a+ = 2g, a. 



(2.40) 



Then the spectrum of allowed charges fl2.39| ) coincides with ( |2.3U[ ) and the dimensions of 
the Kac spectrum ( |2.31| ) read 



(rg - s) 2 - (g - 1) : 

4g 



(2.41) 



The corresponding central charge (|2.27| ) 

C = 1-6(0- lf/g 



[2A2) 



is symmetric under g — > 1/g. 

By (|2.40|) the coupling constant g of the gaussian field is related to the "vac- 
uum" eigenvalue (3 of the connectivity matrix 



(3 = 2cos(7rp ) = — 2cos(7T(7) 



(2.43) 



The branch of the multivalued function g = — arccos(— (3) is determined by the dynamics, 
that is, by the choice of the critical regime for the SOS model. The analysis of the 0(n) 
model []22|1 which is in a sense dual to the SOS model with f3 = n, suggests that the interval 
< g < 1 describes the dense phase and the interval 1 < g < 2 describes the dilute phase 
of the SOS model. Sasha Zamolodchikov conjectured |25| that the other branches of ( |2.43| ) 
( 2 < g < oo) correspond to multicritical regimes of the SOS model. This is shown to 
be the case for the SOS model on a fluctuating lattice ||11||. We will give a sketch of the 
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analysis of ref. [11] in one of the next sections. The m-critical point corresponds to g in 



the interval (m — 1 , m) . 

The order parameters of the SOS model are represented by vertex operators in the 
gaussian theory 



x *(s) = gfcrCp-poM-) <_> V a {Z) = e i " ax ®, a = p- Po (2.44) 



The operators carrying charges p = rpo,r = 0, ±1,... have conformal dimensions in the 
Kac spectrum 



A electric _ V 2 ~ vj _ (g ~ l? (r 2 - I) 

~ 4 g ~ Tg ~ rr 



electric 



The vortex operators in the SOS model can have only integer magnetic charge. The 
conformal dimension of a magnetic operator with charge \i = m,m = 0, ±1, ... is 



A magnetic ( m fi'/2) (g 1) « /- ,„n 

A m = J = A m0 (2-46) 



4g 



The presence of negative dimensions means that the theory is not unitary. There are two 
charge conjugated identity operators (p = ±po). 



2.3. The loop expansion 



The partition functions of the IRF models on a surface with the topology of a sphere 
is equal to the partition function of a gas of nonintersecting loops on the surface [p2[ , . 
The loops have the meaning of domain walls separating domains of constant x. 

The microscopic construction of the loop expansion goes as follows. Write the partition 
function as a sum of monomials, choosing for each plaquette only one of the seven terms on 
the r.h.s. of ( |2.17| ). The ensemble of field configurations S — > X is divided in this way into 
subsets, each characterized by a given decoration of the plaquettes of S. The decorations 
of all plaquettes form a set of nonintersecting polygons on the lattice S* dual to S (Fig. 
3). The polygons play the role of boundaries of the domains of constant x. 
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Fig. 3: A loop configuration formed by decorated plaquettes. 



The partition function is equal to the sum over all configurations of polygons on the 
world sheet and all allowed values of x in the connected domains bounded by polygons. 
The Boltzmann weight of each configuration is equal to T~ L where L is the total length of 
all polygons, times a product of factors associated with the domains bounded by polygons. 

Let us consider for simplicity a lattice with the topology of a disc with and impose a 
Dirichlet boundary condition x = constant. The weight of a domain of constant height x 
with n boundaries is (S x ) 2 ~ n . To prove this we can use the freedom to redistribute the 
Boltzmann weights between plaquettes and vertices. The new plaquette weights are given 
by (glj) with 



X 3 X 4 



X l X 2 



Ox±x s XS X4 0x^X2 0x2X1 



1 



s 



a; 4 



rj-,Cx 1 x 3 Cx 3 X4 ~X4X2°~X2X 1 g 



(2.47) 



2-4 



X l X 2 
X 3 X 4 



n n n n 

^X\ X 3 ^X 3 X4 ^X4X2 ^ X2X\ 



s. 



1 4 



and the new vertex weights are 



Q(x) = S x 



(2.48) 



Using the fact that the number of bonds is twice the number of faces, one can easily show 
that the total power of S x is equal to the Euler characteristic x = 2 — n of the domain. It 
is convenient to distribute the weight S^T n as follows: a factor S x for the outer boundary 
and a factor S 1 " 1 for each of the inner boundaries. 

Next, for a given configuration of polygons we wish to sum over all possible heights of 
the domains bounded by these polygons. Let us start with the domains with the topology 
of a disc, i.e., having no inner boundaries. The sum over the allowed heights of such 
domain is performed using the definition ( |2.2j ). The result is f3S x where x is the height 
of the surrounding domain. The factor S x cancels with the factor S~ l in the Boltzmann 
weight of the surrounding domain, associated with this boundary. Proceeding this way, 
from the inside out, we find that the Boltzmann weight of a polygon configuration is 
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obtained by assigning to each polygon a factor (3 = 2cos(7rpo)- There will be also an 
overall factor S x where x is the height of the most outer domain . Note that the polygons 
have no orientation. One can introduce an orientation and split the weight of a polygon into 
two phase factors exp(±«7rpo) for the two possible orientations. Note that the Boltzmann 
weights of the loops do not depend on the local curvature. 

The correlation function of two order parameters, or vertex operators (|2.9|) is equal, 
up to a normalization, by the partition function of the loop gas where the loops enclosing 
only one of the two points have different fugacity (3 p . The correlation function of two 
desorder parameters, or vortex operators with magnetic charge m is given by the partition 
function of the loop gas in presence of m nonintersecting lines having as extremities the 
two points. We have seen that the partition function on a surface with the topology of a 
disc or sphere depends on the model only through the momentum po of the ground state in 
the target space. The mapping onto the SOS model is trivial in this case. For higher genus 
surfaces the partition function will depend on the spectrum of allowed momenta (torus), 
fusion rules (double torus), etc. The mapping onto the SOS model then becomes more 
and more involved and requires the introduction of a system of distributed electric and 



magnetic charges |26] |T3|]. For more complicated geometries it is convenient to formulate 
the loop expansion using a special diagram technique which will be considered later in this 
section. 



2.4- Summing over the world-sheet geometries 

Now we come to the problem of the evaluation of the string path integral. The 
partition function of the string is defined formally as an integral in the space of all embedded 
surfaces . In our case the integration with respect to the intrinsic geometry is equivalent 
to the average of the partition function in the ensemble of all possible lattices S. The 
dual graphs S* are generated by the perturbative expansion of a 4> 4 matrix field theory. 
In this approach the length of a loop on the graph is an integer (all links of the graph are 
supposed to be of unit length). In what follows we prefer to consider the length L as a 
continuous parameter, in order to simplify the notations. However, the equivalence with 
the discrete formulation is complete. 

For the moment both the length L and the area A are dimensionless quantities. Since 
we are interested in the scaling limit A — > oo, L — > oo , later we are going to introduce the 
renormalized quantities 

£ = aL, A = a 2v A (2.49) 

where a will be a cut-off parameter (elementary length). It is natural to assume that the 
dimension of the world sheet is 2; then 1/v gives the fractal dimension of the loops which 
is determined by the dynamics. This is also the fractal dimension of the boundary of 
the world sheet with Dirichlet boundary condition on x. It can be any number between 1 
and infinity. 

Let S be any closed connected discrete manifold whose cells are squares. Introducing 
the bare "cosmological constant" Kq coupled to the area (= number of plaquettes) A(S) 
and the bare "string interaction constant" kq coupled to the genus (= number of handles) 
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H(S) , we define the canonical partition function as a sum over all such S 



F(K , «o) = 4~ 2H(S) e- KoA ^F(S) (2.50) 



where F(S) is the partition function for frozen geometry defined by ( |2.22 ). This partition 



function is well defined for Kq larger than some critical K* where the continuum limit may 
exist. In what follows the cosmo logical constant will be introduced implicitely through the 
measure in the space of empty surfaces. 

The loop expansion acquires a new significance in the case of fluctuating geometry. It 
helps to reduce the problem to the problem of surfaces without embedding. The condition 
for this is the factorization of the Boltzmann weight of a polygon configuration to a product 
of weights associated with the connected domains. This property is very restrictif. It 
excludes, for example, terms in the action quadratic in the scalar curvature. 

The integration measure over surfaces is controlled by two parameters: the cosmolog- 
ical constant Kq coupled to the area of the world sheet and the string interaction constant 
kq coupled to its topology. If F^ (A) is the partition function of surfaces with fixed area 
A and topology of a sphere with H handles, then the canonical partition function reads 

oo /»oo 

F(K ,k ) =J2 K o dA e~ KoA F^(A) (2.51) 

H=0 J ° 

Note that Kq and kq are dimensionless constants which will be renormalized in the con- 
tinuum limit A — > oo. There are different ways to achieve the continuum limit depending 
on the choice of the measure in the space of surfaces. 

The mapping of the SOS model onto a gaussian field can be constructed also in the 
case of fluctuating geometry. The sum over the geometries means functional integration 
w.r. to the metric g a b in ( |2.23| ). In conformal gauge 



9 ab (0 = gtfiSe*™*™ (2.52) 



where g§ (£) is some fiducial metric, the functional integral leads to a theory of Liouville 



gravity coupled to the matter fields | 2"T|j . As has been demonstrated by David, Distler 
and Kawai , it is consistent to treat the Liouville field <p as a gaussian field (with 
renormalized parameters) and the Liouville interaction as a perturbation. The fields x(£) 
and </>(£) combine into a two-component gaussian field with conformal anomaly 26 defined 
by the action 

A[x, (j>] = - ( d 2 ^detg [irg g%\d a xd h x + d a <j)d h 4>) + £(0(w(0 - £o<t>(0)] 

J (2.53) 
+ A / ^v^et^e 2 ^) 

where A ~ Kq — Kq is the renormalized cosmo logical constant. 
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The vertex operators dressed by the fluctuations of the metric are 



D iir(p-p )x(£)-Tr(e(p)-eo)4>(£,) 



(2.54) 



In particular, the puncture operator ( = identity operator + gravitational dressing) is 
represented by 



no 



r(e(p o )-e o )0(O _ J2tti/0(O 



(2.55) 



The conformal anomalies of the two components of the gaussian field are c x = 1 — 
Ga^/g, Cff) = 1 + Ge\j g and the condition c x + = 26 implies 



4 - Po = 4 f 

We choose the positive solution 

e =g + l, po = g - 1 
so that the two screening charges can be represented as 

a± = po ± e Q = g - 1 ± (g + 1) =$> a+ = 2g, a_ = -2 
The condition that the conformal dimension of the operator ( |2.54|) is one 



(2.56) 



(2.57) 



(2.58) 



+ A, 



P 



Po <P) 2 - 4 



±9 



±9 



combined with ( 2.56 ) leads to the relation 

e(p) 



2 -p 2 =0 



(2.59) 



(2.60) 



which can be interpreted as a mass-shell condition of the same form as the mass-shell 
condition for a light particle in a two-dimensional space-time. 

In principle, the vertex operator ( |2.44| ) dressed by the fluctuations of the metric is a 
linear combination of two operators (|2.54| ) corresponding to the two solutions of ( |2.60| ) 



^(0- A+V p , e (Q + JL.V p ,- e (Q 



(2.61) 



The two terms will have the same dimension if the coefficient A- contains a positive 
power a 2 ' p ' of the cutoff a. Therefore the second term is irrelevant when p ^ 0. Since 
the minimal momentum in the theories with C < 1 is positive (po = 1 all physical 
operators correspond to the positive branch of (|2.60|) 



e(p) = \p\ 



(2.62) 



Another motivation for the choice (|2.62|) is based on the quasiclassical treatment of the 
Liouville theory |28 . 
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In a theory with C = 1 (g = 1, po = 0) both terms become essential in the limit p — > 
and the Liouville interaction in (|2.53|) which comes from the gravitational dressing of the 
identity operator can be not simply an exponential ||29|| . 

By ( |2.57| ) and ( |2.62| ) the Liouville charge of the identity operator ( |2.55| ) in the area 
term of ( |2.53| ) equals Eq — £(po) = 2u where 

_i (s + 1 _ |s _ 1|) = | ?; a»<i; ,, 63) 

The gravitational dimension of the vertex operator with electric charge a rs = p rs + 
Po = r - sg + {g - 1) is ||,fl] 



£o - e{p rs ) 


r — gs 




I<7"1| 


£o - e{p ) 


\9 + l\ 




g-i\ 



In particular, 



(2.65) 



^electric _ M \Po\ ^magnetic _ m ff/2 |go| 



Finally, the so-called string susceptibility exponent 7 s t r is related to the conformal anomaly 
of the matter field by 

C = 1 - 6 7str (2.66) 

1 - 7str 

It can be determined as twice the dimension 5oo in the Kac spectrum of gravitational 
dimensions ( |2.64j ) 



7str = ~2 ; f , ' ~ = -^2Z; v {2 - 7 str) = ^0 (2.67) 





9 ~ 


1 




9 + 1- 


9 


- 1 


V 



2.5. Diagram technique 

We have seen that a map S — > X can be described by a collection of self-avoiding 
nonintersecting loops dividing the world sheet into domains of constant x. The energy 
of each such loop (domain wall) is proportional to its length L. In order to simplify the 
notations, from now we shall consider the length as a continuous quantity. The Boltzmann 
weight of a domain wall is 

T~ L = exp(-2P £) (2.68) 

Further, the weight of each domain of constant height x depends on the geometry of the 
domain through its Euler characteristic. Indeed, the product of the Boltzmann weights 
( p.47|) associated with the sites and squares of the domain is 

(^sites-squares = ^2-2H-n ( 2 _g 9 ) 

where H is the number of handles and n is the number of boundaries of the domain. 

In this way the partition function of an ADE model on a surface S with arbitrary 
geometry has been reformulated as the partition function of a gas of nonintersecting loops 
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on this surface. The Boltzmann weights of the loops depend on the geometry of the surface 
only through homotopic invariants. In particular, they do not feel the local curvature. The 
loops sensible to the topology of the surface are the noncontractible ones that wrap around 
the handles of the surface S. Proceeding as in subsec. 2.3 one can check that all contractible 
loops have fugacity (3 = 2 cos(7rpo) while the fugacity of the loops wrapping around a handle 
is (3 P = 2 cos(7rp) where p is the momentum of the matter-field excitation propagating along 
this handle. Thus the mapping of the ADE models onto the loop gas depends, in the case 
of a surface with arbitrary topology, both on the spectrum of the order parameters and the 
way they interact. This mapping can be systematically formulated using a special diagram 
technique which has been first considered in |16[ and developed further in 0. In the case 
of fluctuating world-sheet geometry this diagram technique gives rise to the Feynman rules 
for the corresponding string field theory, to be discussed below. 

The sum over the embeddings of a world sheet with fixed geometry is a formidable 
problem. It can be avoided by reorganizing the measure over embedded surfaces. First we 
take the sum over the world-sheet geometries with given configuration of domain walls. 
Each domain-wall configuration is determined by its topology and the lengths of the loops. 
It can be described by a Feynman diagram with vertices corresponding to the domains and 
lines corresponding to the domain walls (Fig. 4). 




Fig. 4: A loop configuration on a toroidal surface and the corresponding Feynman diagram 



The heights x and the lengths L are associated with the vertices and lines, corre- 
spondingly. The partition function ( |2.51| ) of the string is equal to the sum of all connected 
diagrams without external lines. 

The weight of a diagram is a product of the weights associated with its vertices and 
propagators. A vertex is defined by the coordinate x, the number of handles H, and the 
lengths Li, ...,L n of the boundaries of the corresponding domain. Its weight factorizes to 
a coordinate- and length-dependent components 

V?°(a;|Li,...,L n ) = {K /S x ) n - 2+2H e- p ^ + - L ^ Wo H \Li, L n ) (2.70) 

where Wo (£i, L n ) is the partition function of nonembedded (empty of loops) surfaces 
with H handles and n boundaries with lengths L±, L n . 
The "propagator" associated with a line L is 

G (x, L\x', L') = C xx > 5(L - L') (2.71) 
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The integration over the lengths is performed with the measure dL/L. This measure is 
consistent with the convention that each boundary has a marked point on it. 

It is convenient to introduce a quantum field theory in which all these graphs appear 
as Feynman diagrams. The latter can be formulated as a functional integral over the 
configurations of a loop field *& X (L) defined in the direct product of the target space X 
and the positive real axis. The interactions of the loop field are given by the vertices ( |2.70| ) . 

It is not difficult to see that the grand canonical partition function Z[J] = exp(i ? [J]) 
(i.e., the one where also disconnected surfaces are allowed) can be written as the following 
functional integral 



Z[J] = J D*e s W 

x,x' 

oo oo -1 poo 

+ EEE^[ dLJ a {L)9 a {L) 



(2.72) 



x n=l H=0 



+ E ^■■ dJ ^vi H \x\L 1 , L n ) n *x(L k 

x J° 1 n k=l 



The multiloop correlation functions are formally obtained as derivatives with respect to 
the source J(L) 



n r 




(9 a , 1 {L 1 )...9 an (L n )) = [ J] 7T7lT F[J] ) (2 ' 73) 



j=o 



The functional integration measure corresponds to the norm 



o 



Formally, eq. ([2.72 ) defines the string field theory corresponding to the target space X. 



However, in order to make this definition more explicit, a lot of work has to be done. First, 
the diagram technique following from the domain wall representation contains tadpoles. 
This means that the string field has a nontrivial background which can be found as the 
solution of the saddle-point equation 

This is the partition function of surfaces with the topology of the disc and bounded by a 
contour of length I situated at the point x. 
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It is easier to solve the classical equation of motion ( 2.75 ) after going to the momentum 
space. For this we expand the string field as a linear combination of eigenstates ( |2.2|) 



Let us denote by Vo (pi, ...,p n \Li, L n ) the vertices ( |2.70| ) in the momentum space. 
They again have a factorized form 

vi H Hpi,..,Pn\Li,..,L n ) = Kr 2+2H N^ Pn e- p ^---+ L ^ Wo II \L ll ... 1 L n ) (2.77) 

where 

n 

^„=En^)(&) 2 — ih 

x k=l 

H 



(p) 



(2.78) 



The p-dependent part of the vertices (|2.78|) has a natural interpretation in terms of string 
states. The factor (S x ) 2 is related to the measure in the X space, each leg is multiplied 
by an order parameter X(p) an d each handle can be considered as the result of contracting 
two legs @ (in the last line of (|2.78|) we have used the orthogonality relations ( |2.11|) ). For 
all target spaces 

N$= 5 (p,p>), N$ pll = Cpp , p „ (2.79) 

where C pp >p» are the fusion rules ( |2.10|) . By applying several times the fusion ( |2.10| ) 
the coefficient ( |2.78| ) can be represented as a sum of products of fusion coefficients. For 
example, 

= Cp p f p f , N^ iP2 = C Pl p 2 pCppf p f , ^p lP2P3 p 4 = C Pl p 2 pCp P:i p 4 (2.80) 
p' p,p' p 

In general, the expression for the coefficient Np 1 /,, Pri can be represented by a "Feynman dia- 
gram" with n — 2-\-2H vertices (to each vertex is assigned a fusion coefficient) and H loops. 



2 Each handle therefore contributes a factor proportional to the volume of the target space. In 
order to define the topological expansion for the string embedded in 2, we have first to introduce 
a cutoff in the dual momentum space. The simplest way to do this is to replace the continuum 
spectrum of momenta by a discrete one with a small spacing 8p. Such a discretization appears 
automatically in the matrix-model realization of the string embedded in 1R [ 30 . The fact that 



the vertices depend explicitely on the cutoff does not contradict the general covariance. Indeed, 
the renormalized string interaction constant k depends explicitely on the cutoff a in the space of 
lenghts I: it vanishes as 1/loga. If we choose 5p ~ 1/loga, then the partition function will be 
cutoff-independent . 
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Thus the p-dependent component of a multiple vertex can be decomposed into elementary 
interactions involving only three strings. Note that all "Feynman diagrams" representing 
the possible decompositions of a multiple vertex have the same contribution (duality prop- 
erty). I 

Now we can write the functional integral ( 2.72j) in the form 



e F ^ = / dVe A W 



dL 

T (2.81) 



i r°° 

■4*] = "2 E / ^( P )(L)(2cos(np))- 1 ^ {p) (L) 
p o 

oo oo _ i poo n 77- 

+ EE E h v<">( Pl , ...,„„!£,, ...Ljn^)^)^ 

n=l fl=0pi,...,p„ ' J ° k=l k 

and the classical equations of motion ( |2.75| ) are equivalent to 

*f p) (L) = 5 P , P0 W(°\L) (2.82) 

x: an poo n it 

By the loop expansion of subsec. 2.3 the quantity (L) is the partition function 

of the gas of nonintersecting loops on a disc with boundary of length L and a fluctuating 
metric. The index (0) means planar topology of the world sheet (H = 0). Eq. ( |2.83|) is 
graphically represented in Fig. 5. 




Fig. 5: A diagrammatic representation of the saddle-point equation. 



Suppose that we have solved the equation for the string backgroung ( |2.82| ). Shifting 
the field in the functional integral ( |2.81| ) by its classical value 



(2.84) 



we arrive at a new diagram technique without tadpoles. The new vertices are obtained 
from the original ones by dressing with tadpoles (Fig. 6). 



We believe that the ^-dependent component of a multiple vertex can be decomposed as well 



into elementary interactions by means of Witten's factorization arguments [31] 
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Fig. 6: Graphical representation of the dressing of the vertices with two and three legs. 



(2.85) 



n -i roc n-\-k 

^)(L 1 ,..., J L n ) = ^i / J] W^(L r )e- 

k=0 ' ^° r=n+l 



PoL. 



e -Po(L 1+ ...+L n ) w (H) {L ^^ L ^ L 



715 lj n-\-l 5 



(2.86) 



Note that for n = 1 ( |2.86| ) coincides with the saddle-point equation ( |2.83| ). The quantity 
W^ H \Li, L n ) is the partition function of the loop gas on a fluctuating surface with H 
handles and n boundaries with the condition that only contractible loops are allowed. 

The gaussian part of the new action (the inverse propagator) is equal to the sum 
of the inverse original propagator ( |2.71| ) and the two-point vertex V^(p,p'\L,L') = 
8p P 'W^(L,L'). (It is convenient to treat the vertices V^ H \pi,p2\L\, L 2 ) with H > 
as part of the interaction.) The propagator G( p )(L, L') is therefore determined by the 
equation 

/AT -i 
— W {0 \L',L)G {p) {L,L") (2.87) 



whose solution is graphically represented in Fig. 7. 



+ <X> + -000+- 

Fig. 7: The full propagator expressed in terms of two-point vertices. 

The symmetric function (L, L') can be considered as a kernel of a symmetric 
operator in the space of the real square-integrable functions defined on the positive real 
line. This operator describes the evolution of the string in the space of the world-sheet 
geometries. Let us introduce its eigenstates which we assume orthonormalized 

/ ^—W ( - 0) {L,L')(L\E) = n(E)(L'\E) (2.88) 
Jo L 
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Then the propagator ( |2.87|) is diagonahzed in the same basis and its eigenstates are 

\ 2cos(7rp) 

The eigenstates (L\E) describe the excitations in the L-space in the same way as the 
eigenstates of the connectivity matrix describe the excitations in the x-space. The wave 
functions of the closed string states are labeled by the "energy-momentum" (E,p) and 
have a factorized form 

(E,p\L,x) = Vfc(E\L) (2.90) 

We have seen that the string background depends on the structure of the X-space 
only through the momentum pg of the ground state (po = 1/h for the Dynkin diagrams 
and po = for the extended Dynkin diagrams) and therefore the same factorization takes 
place for the vertices ( |2.85|) describing interactions due to joining and splitting of closed 
strings. This is the miracle which makes possible to solve exactly the string theory in our 
formalism. Note that in our case the separation of the two global modes occurs before 
taking the continuum limit, which is not the case for the standard discretization of the 
D = 1 string path integral |32|]. It is worth to try to understand better the origin of this 
symmetry. 

Let us mention that there exists a simpler model such that the only allowed momentum 
is pg. This means that the vertices have trivial coordinate part. A microscopic definition 
of this model is given by the 0(n) model on a fluctuating lattice with n = — 2 cos(ng) |33|j . 

2.6. Loop equation 

In order to find the explicit expressions for the dressed vertices ( |2.86| ) we have first 
to solve the saddle-point equation ( |2.82|) . Imagine that the n-loop amplitudes for the 
emply random surface are constructed as the connected correlation functions of an iV x N 
Hermitean random matrix $ with potential U (<&) 

oo n j. 

£ N^ H W^ H \L u ...,L n ) = J] JWl-^J] (2.91) 

H=0 k=l ^ k > 

e F a [J] = J d$e mrC/ ° (,I>) + ^ dLJ(L)tre X p(^) ^.92) 

Then the loop amplitudes Wo H \Li, ...,L n ) satisfy the following system of integral equa- 
tions 



U' {d/dL)wi H) {L,L x ,...,L n ) 

+ E E / dL'W^ H \L';L i \ i£l )W^ H - H \L-L';L j \ j£ j) 

H'=0 I+J={l,...,n} ^° 

n 

W^ H - X \L - L', L ; , Li, L n ) + W { o H) {L + L k ; L s \.^ k ) = 

k=i 



(2.93) 
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Inserting (|2.92j ) into ( 2.83 ) we arrive at the following closed equation for the classical string 
field §, § 

U' o {d/dL)W {0 \L) = j dL'W {0 \L')W {0 \L-L') 

J ° r oo (2-94) 

+ / dL'W^(L')W^(L + L')e- 2PoL ' 
Jo 

which should be completed with the condition that (L) is analytic at the point L = 

OO j- fa 

W^(L) = J2-^W k (2.95) 

fc=0 

In the limit L — > oo the function W(L) is expected to behave as 

W {0) (L) ~ L~ b e PRL (2.96) 

where the entropy per unit length Pr depends on the parameters of the theory. 

Eq. ( |2.94j) has a transparent geometrical meaning []36]. A small deformation of the 
world sheet at the marked point is equivalent to a local deformation of the boundary. Its 
exact form depends on the measure in the space of the world-sheet geometries. In general, 
a local deformation of the boundary is equivalent to a differential operator 

U' a (d/dL) = g 2 d/dL + g 3 d 2 /dL 2 + ... (2.97) 

On the other hand, the deformation of the world sheet cancels part of the integration 
measure in the space of surfaces and produces boundary (or contact) terms. The standard 
contact term, which appears also in pure gravity, is due to the surfaces with degenerated 
world sheet such that another point of the boundary approaches the marked point. Here 
we encounter a new contact term due to one of the loops approaching the marked point 
at the boundary. 

The loop equation Q2.94D has been derived originally fjjin terms of the coefficients W k 



n+k-2 



k>2 



71-2 oo ( , } (2.98) 

fc=0 P,q=0 

We arrive at this form of the loop equation if the gas of loops is considered on a random 
graph, the length of the loops being the number of bonds they occupy. The formulation 
in terms of continuous lengths has been proposed by Kazakov [|36 . 
It is convenient to introduce the Laplace image 

POO 

W(P)= dLW^{L)e- pL (2.99) 
Jo 
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which can be interpreted as the partition function of surfaces with the topology of a disc 
and "boundary cosmological constant" P. In order to simplify the notations we do not 
assign an index (0) to W{P). The condition of analyticity ( |2.95| ) implies 

W{P) = -f + V j± + P^oc (2.100) 

(In what follows we assume the standard normalization Wq = 1.) The asymptotics Q2.96Q 
means that the series (|2.10C1| ) converges for |P| > Pr; in the vicinity of the point Pr it 
behaves as (P - P R ) b ~ 1 . 

We are interested in the solutions of the loop equation having a single cut [Pl,Pr] 
along the real P-axis. The loop equation (|2.94| ) in terms of P-variables reads 



W\P) = $ 2^(P^P') [ ~ U '° {P,) " PW(2P ° ~ P ' )] (2 - 101) 

where the contour C encloses the cut [Pl, Pr] of the Riemann surface of the function W(P). 
The contour integral makes sense only if the two cuts [Pl, Pr] and [2P — Pr, 2P — Pl] 
of the Riemann surface of the integrand function do not overlap, that is, for Po > Pr. The 
positions Pl and Pr of the branchpoints are functions of the potential U (P) and the 



temperature Po of the loop gas. Eqn. (|2.101| ) implies the following conditions on the real 
and imaginary part of W(P) 

2Re#(P) + PW(2P -P) + U' (P) = 0,P e [P L , P R ] 

2.102 

ImW(P) = 0,Pg[P L ,P R ] 



Critical behaviour 



3.1. Critical surfaces, noncritical loops 



The partition function (|2.50|) is defined in certain domain T> in the space of the param- 
eters Po , Kq , . . . . The critical behaviour is achieved at the boundary dV> of this physical 
domain, where the partition function and the observables develop singularities. There are 
two types of singularities in the space of parameters. The first is related to the diverging 
area of the fluctuating surface. It appears when the coupling constants g±,g2, ■■■ in the po- 
tential U (P) which defines the integration measure in the space of nonembedded surfaces, 
are tuned in a special way. The loops on the world sheet are small and do not affect the 
scaling properties. The critical behaviour of the model is the (multicritical) behaviour of 
pure gravity. The m-critical points of the nonembedded random surface |3(| have been in- 
terpreted [j37| , |T38l| , [|39l , PHI as pure gravity coupled to a conformal field theory (2, 2m — 1). 
It can be mapped onto the Coulomb gas ( |2.53| ) with g = m — 1/2. The m-critical potential 
has been chosen in [3^] as an even polynomial of degree 2m 



Uo(p) = J2(- 



fe _ x (fc-l)!m! 
(2fc)!(m — k)\ 



P 



2 k 



(3.1) 
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However, if the symmetry P <-> — P is abandoned, it suffices to take a polynomial of degree 
to + 1. The scaling limit is achieved by replacing L, P by the renormalized quantities £, z 

L = £/a, P = P* + az (3.2) 

At the m-critical point the scaling part w(z) of the planar loop amplitude is defined by 

W{P) = W(P*) + a m - 1/2 w(z) (3.3) 

In the matrix-model realization of the multicritical points the renormalised cosmological 
constant A scales as a~ m which could mean that the fractal dimension of the boundary 
of the world sheet is 2/m. For to > 2 this is less than the classical dimension 1. The 
resolution of this paradox comes from the observation |37|] |38|that for m > 2 A is coupled 



to the operator with minimal (negative) gravitational dimension 8 m i = 1 — m/2 and 
therefore does not measure the area of the world sheet. The true cosmological constant A 
has dimension a~ 2 for m > 2 and a -1 for m = 1. Finally, the string interaction constant 
scales as 

Ko = K a 1 "^/ 2 ; 7str = 3/2 -to (3.4) 



3.2. Dense phase 

When Pq approaches its critical value P* , long-range effects due to the diverging size of 
the loops change the critical behaviour of the model. For generic potential U the random 
surface grows only due to the critical loops which are densely packed on the world sheet. 
Its area is essentially equal to the total length of the loops. This is the dense phase of the 
loop gas. If both the potential U and the temperature Pq are tuned, a more complicated 
picture arises. 

The regime of critical loops coupled to noncritical surfaces describes the dense phase 
of the loop gas model. In this regime the loops are densely packed and fill almost all 
surface of the world sheet. The nonrenormalized area of the world sheet is therefore equal 
to the total (nonrenormalized) length of the loops: A = L to t and P plays the role of a 
cosmological constant. 

It is easy to see that the condition for complete compensation of the entropy and the 
energy of the loops is Pr = Pq ■ First let us remind that W(L) behaves when L is large 
according to (|2.96| ). Let us choose one of the random loops on a surface with the topology 
of a sphere. It splits the sphere into two discs, each contributing a factor e PRL when L 
goes to infinity. Thus the entropy 2PrL of the loop totally compensate its energy 2P L 
when Pr = Pq. 

The vicinity of the critical point is defined by introducing a cutoff parameter a with 
dimension of length. The variable z dual to the renormalized length I (the renormalized 
boundary cosmological constant) is introduced by 

P = P* + az, P R = P*- aM (3.5) 

The parameter M is the contribution of the fluctuations of the surface to the renormalized 
boundary cosmological constant. It defines the position of the cut in the z-plane. 
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The singular part of the loop amplitude W(P) behaves at the critical point Pr = 
P (M = 0) as z 9 where g is the solution of 2 cos(Trfir) + /3 = 0,0<g<l ffi\M- 
Therefore we define the scaling part w(z) of the loop amplitude as 

W(P) = W* + (constant) a 9 w(z) (3.6) 

where W* is the value of the bare loop amplitude for P = P* , Po = P*. 

Now we are going to find the scaling law for the cosmological constant A which is 
defined as 

P = p„ _ (constant) a 2v A (3.7) 

If we assume that the dimension of the world sheet is 2, then the exponent v has the 
meaning of inverse fractal dimension of the boundary of the world sheet. The induced 
boundary cosmological constant M is a function of Po which vanishes at Pq = P*. Therefore 
the renormalized bulk and boundary cosmological constants are related by 

A = (constant) M 2v (3.8) 

The scaling of A can be extracted from a nonlinear algebraic equation which follows 
from the loop equation ( |2.101| ). In order to get rid of the contour integration we use the 
following trick. First write again eq. ( |2.101|) with the variables P, P' replaced by their 



images w.r.t. the reflection P — > 2Pq — P. The integrand remains the same but the contour 
C is replaced by a contour C enclosing the cut [2P — Pr, 2Po — Pl] of W(2P — P). Then 
add the two equations and apply the Cauchy theorem to the integral along the contour 
C + C. The result is the following functional equation for W(P) 

l-U' Q (P)W(P) - P k (9k+2W + g k+3 W 1 + ...) - W 2 (P)} 

k>o (3.9) 
+ [P -> 2P - P] = PW(P)W(2P - P) 

At the point P = Po this functional equation becomes algebraic 

(2 + P)W 2 (P ) + 2U' o {P Q )W{P ) - 2 Po(9k+2 + gk+sWx + ...) = (3.10) 

fc>0 



Eq. (|3 . 1 ) becomes very simple in the limit Pq — » — oo where the only allowed loop 



configurations are those with no space between loops. This extremely dense limit can be 
achieved by choosing gaussian potential (#2 = — 1 5 <73 = ■■■ = 0). Then eq. ( |3.10| ) involves 
only the coefficient Wq = 1 of the expansion ( |2. 100|) and can be solved instantly 



W(P ) = 2/(P + sJP 2 - 2(2 + /?)) (3.11) 

This quantity has a singularity at Po = a/2(2 + ~0) = P* in accord with the exact result 
obtained with more elaborated technique ||41||. 
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By the definitions ( |3.6|) and ( p77|) we find 

W(P ) = W* + (const. )a 9 w(0) = W* - (const.) v 7 ! (3.12) 

By its dimension w(0) ~ M s and eq. ( |3.12| ) implies v = g. ^From now we will fix the 
normalizations of A and M so that 

w(0) = - cos(ng/2)M 9 , A = M 2g (3.13) 

(Note that the sign of w(M) should be negative.) This scaling law persists for generic 
potential U (P), i.e., in a domain D\ of codimension zero in the space of coupling constants 
^2,^3, ••• • Indeed, let us repeat the above argument for a generic potential U (P) = 
g 2 P 2 + (73-P 3 + .... Expanding eq. (|3.10| ) around the point P = P* we find a relation 



between w(0) = — cos(7T(? /2)M 9 and A. In the limit a — > the most singular part of the 
coefficients Wk is the v.e.v. of the puncture operator V = — <9/<9A 

W k ~ a 2 ^ 1 ' 1 ^ {V) ~ {a 2v k) x -i** (3.14) 

Inserting this in (|3.10| ) we find a relation of the type 

(aM) 2g = A k (a 2u A) k + P^A) 1 "^ (3.15) 

k>l 

where Ai : A 2 , ... and B are numerical coefficients. 

In the limit a — > only the smallest power in a will survive. For a generic potential 
all the coefficients A k are nonzero. The strongest singularity is that of the Ai-term and 
the matching of powers gives 

2g = min[2i/, 2p(1 - 7str )] (3.16) 

Sinse the string susceptibility exponent is allways non-positive, this implies v = g. In order 
to fix the string susceptibility exponent, we use the fact that it is related to the dimension 
of the loop amplitude 

K2-7str)-l=£ (3.17) 

Thus in the dense phase 

v = g, 7 s t r = 1 — 1/(7; < g < 1 (dense phase) (3.18) 

In the interval 1/2 < g < 1 the difference P* — Pq vanishes faster than a and the point 
Po in the definition of the scaling variables can be replaced with P*. This is not the case, 
however, if g < 1/2. In this case the cut of the loop amplitude appears at z = — M — a 2g ~ 1 A 
and in the limit a — > the singularity of the loop amplitude is gone to — oo. We have a 
very poor understanding of this case at the moment. Let us only mention that the fugacity 
of the loops is negative and the fractal dimension of the boundary is larger than two. The 
points g = 1/m, m = 1, 2,3, ... have integer 7 str and should describe topological theories. 
A very interesting discussion on the physics at these points has been given by H. Saleur 
g| (see also 0). 

The limiting case g = 1/2 corresponds to the gaussian model. In this case there is no 
singularity at all and the choice of the distance M is arbitrary. 
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3.3. Dilute phase 

Now let us consider the case when both kinds of singularities are present. If the 
potential U {P) is critical (m = 2) , then at the point M = M* both the length of the 
loops and the area of the space between them become infinite. This phase is known as the 
dilute phase of the loop gas. The multicritical potentials will lead, further, to new types 
of critical behaviour of the gas of loops |TT| which we are going to discuss briefly. 



As we have seen above, the dense phase exists in a domain T>\ of codimension zero in 
the space of coupling constants f?i , ^2, — The generic point on the boundary T>2 = dT>\ 
of this domain corresponds to the next critical regime. The simplest way to achieve it is 
to adjust the constant gs in a cubic potential so that the coefficient A\ in (|3.15| ) vanishes. 
Then the condition of matching of powers reads 

2g = min[4z;, 2v(l - 7str )] (3.19) 

If we are in the interval 1 < 1 — 7 s t r < 2, this condition means 2g = 2v{\ — 7 s tr) • Then 
from ( 3.17| ) we find 



v = 1, 7 str = 1 — g; 1 < g < 2 (dilute phase) (3.20) 

Proceeding in the same way , we can achieve the m-critical point by adjusting the cou- 
pling constants #3, ...,g m+ i so that the coefficients A%, A m -\ vanish. The condition of 



matching of powers then gives [11 



2g = min[2mzv, 2u(l - 7 str )] (3.21) 
Combining ( |3.21|) with (|3.17|) we find 



v = 1, 7s tr = 1 — g; m — 1 < g < m (m— critical dilute phase) (3.22) 

Let us assume that the loop gas with given g describes at the critical point a conformal 
field theory coupled to gravity. Then, by the KPZ-DDK argument ||, || its conformal 
anomaly is given by 

l\ (q ~ l) 2 
C = 1 - 6 htT =1-6^ '- (3.23) 

1 - 7str 9 

Thus the interval < g < 00 covers twice the spectrum — 00 < C < 1 of the central charge. 
Note that even if the central charge is symmetric w.r.t. g — > 1/g, the fractal dimension of 
the boundary is not. In the dilute phase it takes the "classical" value 1 and in the dense 
phase it is greater than one. 

The m-critical quantum gravity coupled to a gas of critical loops exhibits a specific 
critical behaviour characterized by the fugacity f3 of the loops. One can say that the 
(2m + 1, 2)-conformal theory of matter coupled to the gas of loops with fugacity (3 produce 
an effective matter field with central charge C = 1 — 6(g — l) 2 / g where g is the branch of 
the function g = —(1/tt) arccos(/3/2) determined by m — 1 < g < m. 
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The change of the critical regime occurs at the integer points g = m, m = 1, 2, 3, .... 
Let us consider the vicinity of the point g = m dividing the m-critival and m + 1-critical 
dilute phases of the loop gas. In this case eq.( |3.17| ) reads 

M 2 = AL4 m (a 2 A) m " 3 + B] x <* (3.24) 

The first term on the r.h.s. vanishes when g < m because the power of a is positive as well 
as when g < m because then A m = 0. However, when g = m both terms survive and our 
previous arguments may not be applicable. The limits g — > m and a — > do not commute 
and one has to solve the loop equation for g = m and then take the continuum limit a — > 0. 



Exact solutions are known for the cases g = 1 |43[ and g = 2 |4J] , to be discussed below. At 
(7 = 2 the the scaling is the same as in the whole dilute phase, A = M 2 . At the point g = 1 
which is the endpoint of the dilute phase, the scaling law receives logarithmic corrections: 
A = (MlogM) 2 . 

3. 4- The loop equation and the string background in the continuum limit 

We are going to solve the loop equation (|2.101|) in the continuum limit. Let us define 
the vicinity of the critical point as follows 

P = P +az 

p R = p Q - a M (3.25) 
W(P) = (2U' o (P + az) - /3U'{P - az))/(4 - (3 2 ) + (const.) a 9 w(z) 

Here we used as a reference point Po instead of P*; we shall see that if 1/2 < gr < 1 this 
makes no difference in the limit a — > 0. 

After the vicinity of the critical point is blown up by the change of variables ( |3.25|) , 
the cut [Pl, Pr] of the Riemann surface of W(P) is replaced by the cut [— oo, —M] along 
the negative z-axis. Eq. ( |2.102| ) implies the following conditions on the real and imaginary 
parts of w(z) 

2Rew(z) + f3w(-z) =0, z < -M 

(3.26) 

Imw(z) = 0, z > -M 



Eq (|3.26| ) becomes more transparent when written in terms of the variable r which we 



have introduced in ||,|| 

z = Mcosh(r) (3.27) 
Then the z-plane cut along the interval — oo < z < — M is mapped into the half-strip 

Rer > 0, -7r < Imr < n (3.28) 

The two sides of the cut are mapped to the horizontal boundaries [r = t ± in, t > 0] of the 
r-strip. Therefore , if z is real and positive, then — z ± iO = Mcosh(r ± in). Since w(z) is 



a real function ( (w(z))* = w(z*)), the first of the eqs. (|3.26| ) can be written as 



[cos(ird/dr) - cos(tt g)]w(z(r)) = (3.29) 
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Eq. ( |3.29| ) has a single solution (up to a constant factor) which is real along the vertical 
boundary of the r-strip and behaves at infinity as z 9 



w(z) = —M 9 cosher) (3.30) 
If we return to the original variable z, the solution reads 

w(z) = -h(z + \/z 2 - M 2 ) 9 + (z- Vz 2 - M 2 ) 9 } (3.31) 

It can be expanded as an infinite series in fractional (in general) negative powers of z with 
radius of convergence 1/z = 1/M 



,( z ) = wtM 9 {M/z) 2n ^ 9 (3.32) 

and dimensionless coefficients 



w( 

n>0 



r(2n±fir) 



T(n + l±g)nl , g ^ 

= ± (2n±g)(2n± g - l)...(n + 2 ± g) 



Even if each of the terms in ( |3.32| ) has a cut — oo < z < 0, the whole series defines a 



function which is analytic for \z\ < M. 

The v.e.v. of the loop operator is determined by the same equation ( [3.26|) both in the 
dense and dilute phases. It is given by the same analytic function ( |3.31| ) with g ranging 
from to oo. The m-critical behaviour of the loop amplitude is described by the branch 
m — 1 < g < m of the parametrization ( |2.43|) . The function w(z) has a square-root 



singularity at z = —M . It is meromorphic in the z-plane with a cut from z = —M to 
z = — oo where it behaves as z 9 . For nonrational values of g the Riemann surface of w(z) 
is infinitely foliated and has two cuts — oo < p < —M and M < p < oo. All sheets except 
the first one have two cuts. We will always consider the function w(z) on the first sheet. 
For g = p/q the Riemann surface of w(z) has a branch point of order q at infinity. Then 
the first and the last sheet have only one cut. 

If we introduce a renormalized length £ = aL, then the inverse Laplace image of w(z) 



ZOO POO 



w(£)= dze lz w(z)] w(z)= d£e- zt w(£) (3.34) 



o 



is the Bessel function [45 



M 9 

w{£) = —K g (M£). (3.35) 



It satisfies the following loop equation 



d£'w{£')w(£ - £') + j d£'w(£')w(£ + £') =0 (3.36) 
(i Jo 
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which is the renormalized version of (|2.94 ). The renormalized string background ( |3.34| ) 



behaves differently at large and small lengths 

^=^M° m ~{f;_T t ~ me - Me > ££Jp (3-37) 



The small-£ asymptotics (|3.37|) can be taken as a subsidiary condition to eq. ( |3.36| ) in 
order to have a unique solution. The large-£ asymptotics describes a loop of length much 
larger than the characteristic length I = 1/M . The area of the world sheet is zero in 
this limit and the power of I is the same as in the topological gravity (g = 1/2). The 
small— £ asymptotics describes the critical point where the area of the world sheet and 
the characteristic length of the loops are infinite. We see that the square-root singularity 
changes generically to a nonrational one in the limit M — > 0. 

In order to to reproduce the regime of critical surfaces, noncritical loops we have to 
take the limit A — > , keeping M finite, i.e., to leave the trajectory A = M 2v ' . However, 
we can achieve the critical points of nonembedded surfaces by tending the fugacity (3 of 
the loops to zero. This limit is achieved at the half-integer values of g. By the scaling of 
the loop amplitude w(z) these points can be identified with the multicritical points of pure 
gravity ||36|1 - The explicit expression for the loop average w(z) in the m-critical point is 

' ^om-l/2 (2m -1)! M \m-l/2-k/2( 9 M \k/2 



w(z) = - V(l + (-) fc )2 w -^ K ^ (z + M) m ~ ' ' (z - M) k ' z (3.38) 

Eq. ( p.38|) reproduces the known results (see, for example, f3"4| ) for g = 3/2 (pure gravity) 
and g — 1/2 (gaussian model) 



-2(2 z-M)7(F +MP, 9 = 3/2 
' -2^/(2 + M)/2, g — 1/2 



At the integer points g = m, m = 1, 2,3,4, ... the planar loop amplitude ( |3.31| ) is a 



polynomial and has no cut at all. However, this is not always the physical solution because 
the limit a — > is taken prior to the limit g — > m. The exact solution for the case g = 2 
4"4j| lead to the following expression in the continuum limit @ 

w(z) g=2 = - ( z ± ^ z2 ~ M2 ) 2 lo g ( z ± V 7 ^ 2 - M 2 

(3.40) 

4-wlz)] , A = M 2 

dg J g=2 



4 The equivalence between the gas of dilute loops with fugacity (3 = —2 (g = 2) and the bosonic 
string embedded in -2 dimensions can be established as follows. If the loops are considered as 
oriented, then their fugacity is -1. Further, it is easy to see that the condition of nonintersection 
can be abandoned, since the total contribution of the configurations containing intersecting loops 
vanishes because of cancellations. Therefore the partition function can be written as the expo- 
nential of minus the entropy of a single oriented loop, with no restriction to its configurations. 
This is exactly the determinant of the Laplace operator in the lattice S, which can be written 
as a Gaussian integral with respect to a couple of grassmanian fields x, x defines on <S. Taking 
the derivative w.r. to g we subtract the zero mode of the Laplacian. By the Kirchoff theorem, 
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The case g = 1 has been solved by M. Gaudin [^3| for the gaussian potential Uq = 
— 1/2<E> 2 . The imaginary part of w(z) along the cut — oo < z < — M has been found as 
the solution of an integral equation with Cauchy kernel. From this it is not difficult to 
reconstruct the meromorphic function w(z) 



w(z) 



-M[t 2 coshr + 2r log(oM) sinhr] 



= -2 



M 2 log(aM) log 



M 



M 2 \ 



log 

d 2 



+ v 7 ? 



M 2 



M 



(3.41) 



w(z) + (const) z, 



3=1 



A = [Mlog(aM)] : 



The origin of the logarithmic corrections in these two cases is quite different. In the case 
g = 1 this is the appearance of a massless excitation (the "tachyon" ) . The divergences due 
to this zero mode produce the logarithmic factor in the relation between A and M. The 
logarithmic scaling violation is a well known phenomenon in the string theory with X = H 

EI- 



it 



In the case g = —2 there are no divergences which can produce logarithmic 
corrections in the scaling of M. The logarithmic factors in the loop amplitude are related 
to the fact that the partition function vanishes and therefore the loop amplitude is actually 
a derivative ||42|| . Within the matrix-model representation, the vanishing if the partition 
function can be explained with the occurence of a supersymmetry at the point g = — 2 
This has nothing to do with the sum over geometries. The partition function 



TT 



for fixed world-sheet geometry is equal to the determinant of the Laplace operator on the 
corresponding lattice and vanishes due to its zero mode. 

As for the higher integer points, g = 3, 4, it is possible that there are no logarithmic 
singularities i These points will be discussed in 

We would like to understand why the loop amplitude of the g = 1 string is equal to the 
second derivative in g of the generic solution. For the moment this is just an experimental 
fact. 



4. The multi-loop amplitudes in the scaling limit 

4-1. An effective matrix model for the dressed vertices 

In order to be able to exploit the functional integral formulation fl2.72f ) of the ADE 



the finite part of the determinant is equal to the sum of connected spanning trees on the world 
sheet. Since each spanning tree defines a dense loop on the world sheet [46], the derivative of the 



partition function can be also expressed in terms of the gas of dense loops with fugacity j3 — > 0. 
We see that the derivatives of the models g = 1/2 and g = 2 are identical and describe the string 
in -2 dimensions. 

5 Saleur [42] considered the topological points g = 1/2, 1/3, ... which are in a dual to the points 



g = 2,3, .... He argued that a self-consistent definition of the theory as a derivative in g exists 
only for g = 1/2. 
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and SOS strings we need explicit expressions for the vertices ( 2.85Q . Using the random 



matrix representation (|2.91|) ,( p.92|) for the bare vertices, the L-dependent part ( |2.86|) of 



the dressed vertices can be calculated by means of a random matrix model with a special 
potential. Namely, the generating function of the dressed vertices 

00 00 ]\j2-2H-n (■ n AT 

W [ J ] = EEV- W(»\L u ...,L n )]Jj(L k )p- (4.1) 

H=On=o U - J fc=l ^ k 

is equal to the free energy of the effective matrix model 

^=/^e~r-~ (4,) 

with non-polinomial potential 

U(<$>) = U ($)+(3 =-e- 2P ° L wM(L)tTe L * (43) 
Jo L 

The potential U(P) is analytic in the the complex plane cut along the line 2Pq — Pr < 
P < 2Po — Pl; it can be represented as an infinite series in $ with radius of convergence 
2Pq — Pr. Its derivative is related to the planar loop amplitude by 

U'(<f>)=f3W i0) (2P -<f>) + U' o (<f>) (4.4) 

Geometrically, the free energy U[J] of the effective matrix model is the partition 
function of the gas of loops on a surface with fluctuating geometry with the condition that 
all noncontractible loops are suppressed. For example, loops going around a handle or a 
boundary are forbidden. 

The one-loop amplitude in the matrix model ( [4.2| ) coincides in the planar limit N — > 00 
with the string background ( |2.82| ) (the partition function of the loop gas on a disc) 

(N-Hre L *) = W {0 \L), N -> 00 (4.5) 



In order to prove this, it will be sufficient to show that ( |4.5| ) satisfies the loop equation 
Q2.94p . Indeed, the planar loop amplitude W(P) in any matrix model is related to the 
potential U(P) by 

ReW(P)=-lu'(P), P L <P<P R 

2 (4.6) 

ImW(P) =0, P < P L or P > P R 

With the choice (fPD of the potential, ( ^6|) is equivalent to ( |2.102| ) and therefore to (|2.94|) . 
We are interested in the continuum limit of the dressed vertices 



^(£ 1 ,...£ n )=a-^ n - 2+2 ^ W( H Xh/a,...,£ n /a), a^O (4.7) 



Below we give a method for their evaluation using the random matrix representation (|4.2|) . 
The idea is to replace the potential ( |4.3|) with a simpler one such that the corresponding 
matrix model will generate directly the scaling limit fl4.7|) of the dressed vertices. 
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Such a potential will represent a function u(z : M) with the following scaling property 

u(z, M) = p a u(pz, pM) (4.8) 



Comparing Q4.6| ) and the functional equations ( ^.26| ) for the renormalized classical 
loop field w(z) , we see that a possible choice of the scaling potential ((Of) is 



u(z) = (3 



d£ 

T 



w{£)e z 



d_ 

dz 



u(z) = (3w(—z) 



(4.9) 



Let us mention that the inverse problem of finding the potential given the loop amplitude 
has many solutions. For example, the multicritical points of pure gravity can be obtained 



from the polynomial potentials of Kazakov |36| as well as from the non-polynomial singular 
potentials of Gross and Migdal (49| . 

We have shown that the dressed vertices in the scaling limit ( f4.7| ) can be calculated 
as the multi-loop correlators in the effective matrix model with potential 



oo „ n 

W^ H \£ U ...J n )N 2 - n - 2H = / d ^j °°(di/i) W (e)treMm JJ tre 4* (410 ) 
H=0 k=l 



The one- and two-loop correlators in the planar limit are 



w {0) (£) = -M 9 K g (M£) 



1 



(4.11) 



i + u 



-M(l+l') 



(4.12) 



The general formula for the planar multi-loop correlators in the matrix model [50| [Sol reads 



w (0) (4,...,C) 



(£! + .. J n ) V dy 

Qn—2 poo 



2 e -f(y;M)(£ 1 +...+i n ) 



dy 



n-2 



dze - Z (£ 1+ ...+ £n ) 



(4.13) 



f(y;M) 



where the derivatives are taken at the point y = and the function f(y; M) is determined 
by the diffusion equation (Appendix A) corresponding to the potential (|4.9|) 



y= lim (3 / d£ w(£)(y\e z{ - nt)+(d/dt)2) \y) 



d£ 

T 



(3M 9 K g (M£) 



-f(y;M)£ 1 
~7T 



(4.14) 



We assume that the infinite constant produced by the singularity at small distances {£ ~ a) 
is subtracted so that 

/(0;M) = M (4.15) 
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The finite part of the integrai can be extracted by means of analytic continuation w.r. 
to a regulating parameter a. First we multiply the integrand by £ Q+1 / 2 , then calculate the 
integral ( |4.14| ) in the domain Rea > g where it converges, using the formula |5lJ 



POO 

/ di r- x e- zi K a {Mt) 
Jo 



n fur , ^M/2-a r ( Qi + g) r ( Q -g) 



:(M + Z) 



2M r(a + l/2) 

M — z {4A ^ 
2^(1/2 + g, 1/2 - g- 1/2 + a; -^-) 

~, ur , M/2-a ^ F (V2 + S + fc)r(l/2 — g + k) / M — z^ k 



2M { ' ^ T(l/2 + a + k)k\ V 2M 

and finally perform analytic continuation to a = 1/2. 
The result is 



y=/3M V2M (M + /) ^ 3H(3fc3i)! 

=2(27r) 3 / 2 (^ 2 - 1)M« +1 / 2 (4.17) 

1 r(fc — 1/2 + g)T(k — 1/2 — g) / M — f 
f^kUk-Dl T(1/2 + q)T(1/2-q) V 2M 



(the = term disappears because of the T(0) in the denominator.) 

The parameter y is the coupling constant of the matrix model; in terms of the loop 
gas it is the coupling of a special local operator, to be discussed below. 

4-2. Some critical exponents 

Until this moment we have considered two of the critical exponents of the loop gas 
model - the string susceptibility 7 s t r and the fractal dimension of the loops Db = 1/v 
which is the same as the fractal dimension of the boundary of the world sheet. Now we 
are going to introduce a third exponent Dc which has the meaning of fractal dimension 
of the connected domains bounded by loops. 

By construction, the vacuum energy of the matrix model in the planar limit is the 
partition function of the gas of nonintersecting loops on the sphere, with one labeled do- 
main. One can imagine it as the partition function of the gas of small spots (baby universes) 
floating in the labeled domain. A spot of perimeter I represents the partition function of 
the loop gas on a disc and have a Boltzmann weight u(£) = f3w(£). The size of these 
baby universes is of the order of the cut-off a because of their entropy being divergent at 
I — » 0. The measure in the space of world-sheet metrics is tuned to the m-critical phase of 
nonembedded random surfaces. The parameter y makes sense of "cosmological constant" 
for the labeled domain. 
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Let Ac and A be the area of the labeled connected domain and the total area of the 
world sheet, correspondingly. (Of course, the labeled domain is equivalent to any other 
connected domain on the world sheet.) Then the fractal dimension of a single connected 
domain is determined by the way these two areas grow near the critical point M = 

A c ~ A Dc ' 2 (4.18) 

The area of a single connected domain is measured by the matrix- mo del puncture operator 

V c = -j- (4.19) 

The coupling constant y of the matrix model scales as M 3+1 / 2 ; this is sufficient to determine 
the fractal dimension Dc- It is convenient to consider the dense and the dilute phases 
separately. 
(i) Dense phase 
In this phase 

Va = -^A C ,V=^A (4.20) 

The cosmological constant A is related to y by y = M 9+1 / 2 = A 1 / 2 " 1-1 /^ 3 ) and the fractal 
dimension of the connected area is 

Dc = l + 1 (4.21) 

For positive (3 (g > 1/2) all Boltzmann weights are also positive and the gas of loops 
allows a statistical interpretation. Then the area of a connected domain diverges more 
slowly than the total area 

A conn ~ (4.22) 

In the dense phase A coincides with the total length of the loops on the world sheet, 
and A conn - with the total length of the loops forming the boundary (hull) of a connected 
domain. 

The scaling ( |4.22| ) is confirmed by the following argument based on the Coulomb gas 
picture. 

The susceptibility d 2 U/dy 2 of the matrix model can be interpreted as the two-point 
function of a special local operator defined on the whole world sheet and not only on the 
labeled connected domain. The two-point function of this operator is zero when the two 
points are in different domains and one if they are in the same domain. This operator 



can be constructed as a vertex operator with electric charge e = g — 1 ± 1/2 ||52|| . The 
corresponding gravitational scaling dimension 

<W/2 = 1/2 ~ 19 ~ M = 1/2 - 1/(40) (4.23) 
2g 

is positive when 1/2 < g < 1. The dimension 1 — 5o,i/2 = Dc/2 of the coupling constant 
of this operator indeed coincides with the dimension of y. 
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(ii) Dilute phase 

In this phase the puncture operator Vc of the matrix model scales as A~( ff+1 / 2 )/ 2 . In 
the m-critical phase Vc = A™^ 2 and the two areas are related by 

A c ~A Dc/2 ; D c = (2g + l)/m (4.24) 

In the Coulomb gas picture, the constant y is coupled to the same local operator with 
gravitational scaling dimension 

*o,i/2 = 1/2 ~^~ 11 = 3/4 - g/2 (4.25) 

This time 1 - <5 ,i/2 = (m/4)D c . 

Knowing the area of a connected domain, we can evaluate the characteristic number 
of connected domains on the world sheet. By the Euler theorem, this is also the number 
of loops N (including the most external loop representing the boundary). Therefore 

Af = A/ A c ~ A l ~ Dc 

A i/2-i/(4 5)) denge phage (4.26) 

j^(m- g -i/2)/m^ m-critical dilute phase 

In the interval 1/2 < g < 3/2 where the loop gas has statistical interpretation (all Boltz- 
mann weights are positive, including these related to the measure over random surfaces) 
the total number of loops tends to infinity when the critical point M = is approached. 
The density of loops is proportional to 1/Ac and goes down to zero at the critical point. 

The number of loops becomes one at the half-integer points g = m— 1/2, m = 1, 2, 3, ... 
which are the solutions of the equation (3 = 0. These are the multicritical points of pure 



gravity |36| where the only loop is the boundary of the world sheet. At m = 1 and m = 2 
the dimension <5 0)1 / 2 of the local operator coupled to y is zero and therefore Vc = V. In the 
higher multicritical points with (3 = the dimension of the operator Pc fits the dimension 
5^-1,1 of the Kac table of the model (2, 2m — 1) 

l/2-|m-3/2| (m-l)-(m-l/2)-|m-3/2| 
°o,i/2 = 2 = 2 = m " 1 ' 1 ( 

Therefore Pc can be identified with the operator with the most negative dimension in the 
m-critical regime of pure gravity [JJ71 |OT . 



4-3. Feynman rules 

Now we are able to formulate the Feynman diagram technique for the ADE strings 
in a more explicit way. 

The vertices (|2.77Q and (|2.85[) read, in the continuum limit 



v^ H \x\£ 1 ,...J n ) = (kS x ) 2H 2+n w^(£ u ...,£ n ) (4.28) 



38 



u^Cpi, -C) = ^-^iV^^W^, ...,* n ) (4.29) 

where are the loop amplitudes in the one-matrix model with potential u(z) = (3w(—z), 
w(z) given by ( |3.31| ), and the renormalized string interaction constant k is defined by 

k = a- ui2 -^K = a-^Ko (4.30) 

The explicit form of the vertices without handles (H = 0) is given by ( |4.17| ) and Q4.13|) . 
To our knowledge, the explicit formula for the general macroscopic loop amplitudes in the 
one-matrix model is not yet found. Below we give an argument allowing to imagine the 
general form of the matrix-model loop amplitudes with H > 0. 

In the general case the loop amplitudes in the effective matrix model can be con- 
sidered as loop amplitudes for the gaussian model dressed by tadpoles representing the 
nongaussian part of the potential (see eq. ( |2.86| )). From the explicit form of the gaussian 
loop amplitudes (Appendix B) we conclude that 



w^^.^Q = V / U;Pf ) (4 + .Jn)e- M(U " +4) (4.31) 

where PA {£) is a polynomial of degree n — 3 + 3H with coefficients depending on g and 
M 

n-3+3H 

PL H \Z)= E Ai "k(9,M)t k (4.32) 

k=0 

(Note that n — 3 + 3H is the complex dimension of the moduli space for surfaces with H 
handles and n punctures). 

Let us denote the n-loop amplitude in the continuum limit by 



k=i 



where ^ X {1) is the loop operator creating a boundary of length I at the point x. For the 
tree approximation (H = 0) we will use the symbol (...)o- The v.e.v. of the loop operator 
is, in the tree approximation, 

(*,(£)) = jM9K g {Ml)S x (4.34) 
The two-loop correlator satisfies the equation 

/"OO 

(y xi (£ 1 )<S> X2 (£ 2 )} = Y, C *ix / ^» (0) (4^*x(^ l8 (4)) (4.35) 

Je=o 

where w/°) (£,£') is the two-loop amplitude in the effective matrix model. The latter is 
diagonalized by Bessel functions |)T 



wM = ^e"«<«') = jy Em -L- {m (4.36) 
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(£\E) = -WirE smh(irE)K iE (M£); K tE (M£) = dre- IVUcostlT cos(Et) (4.37) 
7T Jo 



oo 

Mi cosh 



where the states (£\E) form a complete orthonormal system On the other hand, the 
connectivity matrix is diagonalized by 

C xx . = X)Vg ) 2oos(7rp)V^ ) (4.38) 
p 

where the sum includes the allowed momenta p = m/h. Thus the wave functions of 
the closed string states diagonalizing the two-loop correlator are labeled by the "energy- 
momentum" (E, p) and have a factorized form 



(E, p\£, x) = V ( x p) y/E wnh(nE)K iE {M£) (4.39) 
The equation for the two-loop correlator becomes algebraic after the diagonalization 



and its solution is [[10 1 



f°° 1 

For p — 1/2 (pure gravity) and p = 1/3 (Ising model) the two-loop correlator has been cal- 
culated from matrix models (45| , |53| . One can check that (|4.40|) reproduces the expressions 
found for these particular cases. 

Remarkably, the two-loop correlator is universal; it depends on the model only through 
its spectrum of allowed momenta. This makes possible to map the Hilbert spaces of the 
ADE strings onto the Hilbert space of the SOS string whose spectrum of excitations covers 
the whole interval — 1 < p < 1. The two- loop amplitude is essentially the propagator in 
the string diagram technique. The latter is given by the continuum limit of (|2.89f) 

/>oo 

G (p) (£J')= / dE(£\E)G(E,p)(E\£') (4.41) 
nt-p \ 2 cos HkE) cos(ttp) 

G{E,p) = — — — r — r (4.42) 

cosh{nE) — cos{7ip) 

The (E : p) space has the geometry of a semi-infinite cylinder. It is represented by the 
infinite (E : p) plane factored by the relation of equivalence 

(E,p) = (E,p + 2) (4.43) 

The periodicity of the propagator ( |4.39|) is a direct consequence of the discreteness of the 
coordinate space X. 

It follows from the explicit form of the eigenstates (|4.37 ) that the r-parametrization 



of the boundary cosmological constant z (eq. (|3.27|) ) provides a local coordinate dual to 
the energy E. 
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The configuration space of the string fieid theory is a direct product of the X and 
r spaces. The relation between the I and r representations has been discussed in ref. 
||53|| . It has been argued that the r variable makes sense of the time variable in the Das- 
Jevicki approach, and that the transformation from I to r is analogous to the Backlund 
transformation in the Liouville theory. This transformation reads explicitely [|53 



/■oo 

K iE (M£)= dre- MicoshT cos(Et) 
Jo 

f°° rlP 

/ ^ye- MicoshT K xE (M£) 
Jo " 



cos et r°° di ( 4A4 > 



EsmhnE 

and does not relate delta-function normalized bases of wavefunctions. 

Let us consider the case when the target space is the discretized real line 2 C IR. The 
inverse propagator ( |4.42|) acts in the (x, r) space as a differential operator of infinite order 

G' 1 = [coshind/dx)}- 1 - [cosid/dr)}- 1 (4.45) 

This operator can be also interpreted as a finite- difference operator in the x and it di- 
rections. Thus Minkowski rotation r — > t = it of the TL string is described by a lattice 
Hamiltonian. The saddle-point equation ( |3.29|) means that the classical string background 
is annihilated by the Hamiltonian ( |4.45| ). 

It will be convenient to apply the Feynman rules directly in the (E,p) space. The 
Fourier image of the ^-dependent part of the vertices ( f4.31|) can be found using the inte- 
gration formula 

Jo t lE{ ] (2M) fc+1/2 r(* + i) {AAb) 

which follows from ( 4.16 ). The result is 

rt H \p 1 ,...,p n \E 1 ,...,E n ) = H J n - 2+2H N Pl ... Pn w i ^ H \E u ...,E n ) (4.47) 



Jo H 4 



' k \e k ) w w{£ u ...,e n ) 



n-3+3H 

= E E m\AW(g,M) (4.48) 

m=0 777,i + . ..+m„=m; 77ij.>0 

fr r- . - , - , r(l/2 +m k +iE)T(l/2 + mk -iE) 

x y/E k smh {i\E k ) f— ^ 

£Ji m k \m k \ 

The vertices in the (E,p) space are given by the product of the r.h.s. of ( |4.48| ) 
and ( |2.78| ). Each vertex represents a symmetric polynomial in E±, ...,E n times a factor 
a/ nEk smh(nEk) / cosh(TcEk) for each line. 



41 



It is convenient to absorbe the non-polynomial factors in the propagators; then the 
new propagator and vertices have the following form 



^ A/sinh(7r£') 2cos(np) cosh(nE) y< smh(n E) /E 

VnE cosh (txE) cosh (txE) - cos (up) cosh(7ri?) 

= Gv{E,p)-Gv{E,l/2) (4.49) 

„ ,„ v sinh(7rS) 2 >^ i , A - n \ 

Gv{E,p) = K -— = }2 T?2^( ^2 4 - 50 

ttE cosh(7ri?) — cos(7rp) n r^ X) E z + (p + 2n) z 

n-3+3H 

vJ£2(pi, ...Pni^i, ..,E n ) = N^ Pn J2 E 

m=0 mi + ...+m„=m; mj,>0 

n rp m k -l 

ml 45 II T% II \ E l + ( s + V2) 2 ] (4.51) 



With the Feynman rules ( fl.49| ), ( |4 . 5 1|) the evaluation of the loop amplitudes can be done 
as in an ordinary quantum field theory i. The Greens functions of this effective field theory 
are the correlation functions of the Fourier-transformed loop field 

/°° rlf 
=X( P ,«M«) (4-52) 



*(B,p) m = VE .^*V *{B,p) (4.53) 
y smh(itE) 

Let us calculate, for example, the three-loop correlator of the loop field ty( p \(£) = 
Y^ x X?p)^x{£)- We start with the corresponding Greens function in the (E,p) space which 
is given by a single Feynman diagram 

3 3 

(J] V(Ek,Pk))o = C PlP2Pa M~^ J] G(E k:Pk ) (4.54) 

k=l k=l 

Then we integrate w.r. to the E variables associated with the three legs according to the 
formula 

Esmh(nE) 
cosh(7r£') — cos(7rp) 
to find 

3 3 

(H <f (Pfe) (4)> = A-(i_w2)c wp 2 P3 n ^i-iai ( m ^) ( 45e ) 

fc=l k=l 

In the particular case Pi = P2 = P3 = Po corresponding to three loops associated with the 
identity operator, eq. ( |4.56| ) is in accordance with the the formula conjectured in 0. 



dE —^j-^ v J t _^ K lE {Mt) = MlK x _\ p \[M€) (4.55) 



6 The form of the propagator ( 4.49 ) suggests an interpretation of the string field as a collection 
of two modes: a bosonic relativistic particle with compactified momentum space in the space-like 
direction and a ghost excitation (twist) with spectrum of momenta p = ±1/2 
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4-4- Spectroscopy 

The spectrum of on-shell states is given by the poles of G(E,p). They form the light 
cone in the (E,p) cylinder (Fig. 8) 

E = ±ie n (p), e n (p) = \p + 2n|; n = 0, ±1, ±2, ... (4.57) 

The periodicity in the momentum space is a consequence of the discreteness of the 
coordinate X space. If the X-space is compact, then the momentum space is discrete 
as well and the on-shell states form an infinite two-dimensional lattice on the space-time 
cylinder. 




-1 1 

Fig. 8: The light cone for a compact momentum space. 
The edges of the strip correspond to identical momenta p = ±1. 

Each on-shell state creates a "microscopic loop" on the world sheet and corresponds 
to a local scaling operator. The sequence of poles corresponding to a given momentum p 
defines the set of local operators 0^ p j. n ; n = 0, ±1, ±2, ... having nonvanishing correlations 
with the vertex operator 



where V x is the operator creating a puncture in the world sheet with coordinate x. These 
operators are the coefficients in the expansion of the off-shell loop field ^^(£) as a series 



of on-shell wave functions [45 



E (\p\+2n)M\W+ 2n \l M+2nl (M£)0 { 



p); 



(4.59) 



n= — oo 



The expansion of the two-loop correlator as a sum of products of on-shell wave func- 
tions reads 

00 II o 

(*(p)(i)*(p)(0>= E -^^ K \ P \+2n(M£) I^iMi'), £'<£ (4.60) 



n= — oo 
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Eqs. ( p39[) and (|Q0|) imply 



{0 {p) . n V {p>) {l)) = M^ +2n \K lpl+2n (M£) (4.61) 

A direct derivation of (Laplace transform of ) ( f4.61|) using the loop equations can be found 
in i and!. 

The gravitational dimension 5 n (p) of the local operator V( p ) ;n can be found from Q4.61] ) 

M \\ P \+2n\ = A S n ( P )-^/2 ^ 5M = \\p\+M-\p Q \ (4 _ g2) 

where po = g — 1 is the background momentum (|2.57|) (see ( |3.18|) , (|3.22|) ). 

The dimensions of the operators Oj^ p y n with p = kpo , k - integer, are contained in 
the KPZ-DDK spectrum ( ^Bl] ) 

S(O rpo;n ) = 5 k +2n,k (4.63) 



4-5. Partition function of the noninteracting string 

The Feynman rules in the continuum limit are generated by the functional integral 

e F[J] = J dq)e A W 

1 f*OC jg rOO -i -i . 



(4.64) 



n=2 H=0 pi ,...,p 



pi H) (^i+...+on*(p h) (v^ A - 



A- J v ^ fc 



which formally defines the string field theory with target space X. The partition function 
of the free string (topology of a torus) is closely related to the vacuum energy of this 
field theory. 

The loop expansion for the embedded surfaces with the topology of a torus involves 
loop configurations with contractible and noncontractible loops. It is easy to see that the 
logarithm of the determinant of the kernel G xx >(£ : £') is equal to twice the contribution of 
the surfaces in the loop expansion containing at least one noncontractible loop. Indeed, if 
we return to the original diagram technique (subsec. 2.5) and retain only the gaussian part 
of the action, the contribution of the gaussian fluctuations of the \l/-field to the vacuum 
energy will be given by the sum of all Feynman diagrams with the topology of a tree with 
one cycle. Each cycle corresponds to a closed path in X and consists of even number of 
points (bare vertices). A cycle with 2k points should be taken with a symmetry factor 2k. 



44 



Let us denote by F\ k . p ] the contribution of the toroidal surfaces (H = 1) containing 
exactly k noncontractible loops with momentum p. Then the partition function for fixed 
momentum p is equal to 

oo 

F (p) = F IM 

The contribution of the surfaces without noncontractible loops is equal to the trace of 
the identity operator in the X-space (= the number of its points) times the l/iV 2 -correction 
to the free energy of the effective matrix model considered in the beginning of this section. 
A simple calculation (see, for example Appendix C of [[49|| ) leads to 

= 9 —^T- log(aM) (466) 

The integral over the length I in (|4.65|) yields a factor log[l/(aM)] (this is the di- 
agonal value of the kernel of the regularized identity operator in the -E-space) and the 
-E-integration gives 



oo 2 







COSh \Kti)' £ 



Inserting this in ( |4.65| ) we find for contribution to the partition function of an excitation 
with momentum p 

F(r) = w »-'^y-M)' ( , 88) 

Let us give some examples. 

(i) The SOS string compactified on TL^h 

Summing over the allowed momenta p = m/h; m = 0, ±1, ±(h — l),h in the 
embedding space X = Z 2 /j, we find 

FM'J) = log(aM) (^±? + 2h i^l) = .og(aM) (»*^*) (4.69) 

The ~E-2h string is well defined if the background momentum po belongs to the spectrum 
of allowed momenta. Therefore the possible values of g are 

k 

g = m ± — ; m = integer, < k < h (4.70) 
h 
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The partition function of the string embedded in ~H-2h is symmetric under the duality 
transformation 

g -> 1/g, h -> 1/h, M — > M 9 (4.71) 

relating the dense and dilute phases. Due to the different scaling in these two phases, the 
cosmological constant is invariant under the duality transformation. The self-dual point is 
g = 1, h = 1. It corresponds to the level-one representation of an SU 2 current algebra. The 
Kosterlitz-Thouless transition may occur at g = l,h = 2 and its dual point g = 1, h = 1/2. 
The case h = 2 describes the model which is a particular case of the Ashkin- Teller 
model. At this point the magnetic (vortex) operator with discontinuity m = 4 (= the 
period in the X space) becomes marginal: its conformal dimension (|2.46| ) is 1. The dual 
model is formally the 7L\ model having as a target space a graph with one point and one 
loop. The space of paths in this space is not trivial because each step is made by traversing 
the loop in one of the two possible directions. This 7L\ string is equivalent to the n = 2 



limit of the O(n) string introduced in fl33|. The vortex operators in the 0(n) string are the 
electric operators in the 7L\ string. The spectrum of momenta consists of all integers and 
the marginal operator is the vertex operator with p = 2 (= the period in the momentum 
space) . 

(ii) The RSOS strings 

The embedding space X = A^-i has spectrum of momenta p = m/h; m = 1, 2, h—1 
and the partition function reads 



(h-l)(h-2) , , 
h ~ 1 12h lQg(aM) 

/# ~ 2 log(o 2(1 - 1 /^A), dense phase (g = 1 - 1/h) 



24 

h(h-l), , 
F Ah ^ = \ 2h M og(aM) 

log(a 2 A), dilute phase (g = 1 + 1/h) 



(4.72) 



(4.73) 



24 

In both critical regimes the partition function of the A^,_i string is equal to the difference 
between the partition functions of the TLiu and 7L<i strings 

F Ah _ t = \{F 2h {g = 1 ± 1/h) - F 2 (g = 1 ± 1/h)) (4.74) 



The relation ( 4.74 ) takes place even before performing the sum over the world sheet ge- 
ometries 0. 

In a similar way one can calculate the partition functions of the D and E strings. 
In the A and D cases the expressions will coincide with these obtained from the matrix 
models in the formulation of M. Douglas by Di Francesco and Kutasov ||55|| . All these 
partition functions can be expressed as linear combinations of partition functions of Z n 



strings 0. The formulas are the same as these for the regular lattice ||12|| . The Coulomb 
gas calculation of the genus-one partition function of the ADE strings was presented in 



Here we restrict ourselves to the unitary theories with < g < 2. 
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5. Concluding remarks 

We have shown that the propagation and interactions of strings embedded in a discrete 
one-dimensional target space are described by string field theory with propagator ( |4.42| ) 
and vertices (|4.48| ). We were not able to calculate explicitely the vertices with higher 
topology but they are defined unambiguously by the matrix integral Q4.2|) . 

Our diagram technique based on the loop-gas representation of the IRF models, pro- 
vides a natural decomposition of the moduly space into elementary cells. The vertices 
represent the topology-changing amplitudes involving no propagation in the target space. 
The fact that the vertices are themselves loop amplitudes guarantees the "stringy" large 
order behaviour [ 57| . The vertices included in a Feynman diagram cut a number of holes in 



the moduli space such that the rest of it is a direct product of one-dimensional spaces (one 
for each propagator). The degenerated surfaces appear when the proper time for some 
of the propagators becomes infinite (the proper time is measured by the number of the 
noncontractible loops along the corresponding cylindric surfaces). This decomposition of 
the moduli space is well defined only if all proper times are nonzero (i.e., if there is at least 
one noncontractible loop in each of the channels. The situation when one of the proper 
times vanish is taken into account by replacing the two involved vertices with a new vertex 
of higher topology. Thus the interaction of the closed string is described by a collection 
infinitely many vertices, one for each topology. This is possibly a general feature of the 
perturbative expansion for any string theory |58[ . 



It would be interesting to compare our diagram technique for the loop amplitudes to 
the one for the string embedded in IR. The two-loop correlator for the string embedded in 
IR il[]5|], n. 

G(E,p)= f 1 2 (5.1) 

Sinh(7T£/J £/ + p A 

is compatible with Q4.40| ) only in the limit of small p and E. The propagator of the 
IR-string contains an infinite sequence of poles at integer imaginary values of E which have 
been interpreted as the special stationary states i discovered by A. Polyakov [K0j. These 
states are created by infinitesimal motion in the X space; in the Coulomb gas picture they 
are constructed by taking the derivatives of the x-field. Therefore, such states should not 
exist in a string theory with discrete target space. The special states in the IR-string are 
related to the closed classical trajectories in the imaginary time direction. 

The infinite spectrum of states in our case comes from the fact that the discreteness of 
the X-space leads to periodicity in the momentum space. As a consequence, the light cone 
maps an infinite set of energies to the same momentum. These energies make the tower 
of gravitational descendents of the vertex operator with given momentum. The classical 
motion in the imaginary time direction is gouverned by a finite-difference Hamiltonian. 

The fact that the difference between the embedding spaces 2 and IR survives in the 
continuum limit can be explained. If we construct the D=l string as an infinite chain 
of coupled matrices |Jl| j62[], then the IR-string will describe the low temperature phase 



8 Since these "states" can be elliminated by redefinition of the vertices, their physical meaning 
is not very clear to us. 
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of the matrix chain and the Z-string will describe the Kosterlitz-Thouless point where 
the distance between two nearest neighbours on the chain remains finite in the continuum 
limit. The change in the critical behaviour of the matrix chain near the Kosterlitz-Thouless 
point is due to the liberation of the angular excitations |63| [B4| 
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Appendix A. The method of orthogonal polynomials 

The integral over Hermitean N x N matrices 

e N2p = [ d$e NtrU W (A.l) 



can be written as an integral w.r.t. the eigenvalues n , n = 1, N, of the matrix $ 

N 



e N F 



k=l i<k 



This integral can be interpreted as the scalar product of the wave function of N fermions 
with itself. The one-fermion wave functions are of the form 

{<j>\n) = P n [<j>)e u W 2 (A.3) 

where P n ((p) ls a polynomial of degree n. The polynomials are fixed by the condition of 
orthonormality 

(m\n) = [ d^Pm^Pn^e 11 ^ = 5 mn (A.4) 



The ground state of the system of fermions is equal to the antisymmetrized product of the 
states \n),n = 1, N 

S(fc,...,tor) = E (-iWl^.-W™) (A.5) 

ctESn 

The operator $ is represented in the space of one-fermion states is represented by the 
tridiagonal Jacobi matrix which we denote by the same letter 

$|n) = R n+1 \n + 1) + R n \n - 1) + S n \n) (A.6) 

The nonzero matrix elements S n and R n = A n _i/A n of this matrix are fixed by the 
recurrence relations 

(n\<Z>U'($)\n) = (2n + 1)/N (A.7) 
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(n\U'($)\n) = 



(A.8) 



which can be interpreted in terms of diffusion of a particle with n-dependent hopping 
parameter. In the limit N — > oo we can replace R n , S n by 



R = lim R n , S = lim S r , 



(A.9) 



and the resolvent of the Jacobi matrix is given by the propagator of a one-dimensional 
relativistic particle 



K nn ,{P) = (n\ 



1 



\n) 



P — S — y?(P — S) 2 — AR 2 



(A.10) 



2R 



v/(P - S) 2 - AR 2 



Another useful quantity is the inverse Laplace transform of (|A.10| ) 

\e L *\n') = e SL I\ n _ n ,\{2RL) 



{n\ 



(A.11) 



where I n (z) is a modified Bessel function. All quantities characterizing the radial part of 
the random matrix can be expressed in terms of QA..10 ) or (|A.11| ). For example, the one- 
and two-loop correlators are given by 



N N 

W(L) = - Y,(n\e L "\n) = e SL - £ I (2R n L) 



n=l 



n=l 



W(L,L')= (n\e L *\ri)(n'\e L '*\n) 

n>N,n' <N 

In the limit iV — > oo we can introduce the continuous variable 

v = 1 

y N 

and consider R n and S n as functions of y: 

R n = R{y), S n = S(y); R(0) = R, S(0) = S 



(A.12) 
(A.13) 



(A.14) 



(A.15) 



These functions are determined by ( |A.7| ) , (|A.8[) . Using the explicit expression ( |A.10| ) for 
the resolvent we find 

■Pn ^ pu , {p) 

< Pl ^{P-S) 2 -AR 2 



tt(1 - y) 



Pr 



dP 



U'(P) 



y/(P - S) 2 -AR 2 
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(A.16) 
(A.17) 



The positions of the branchpoints of W are therefore given by 

P L = S-2R, P R = S + 2R (A.18) 

It follows by (T02|) that 

W(L) = - I dy e s ^ L I (2R(y)L) (A.19) 
Jo 



W(P) = - dy , (A.20) 
Jo ^{P-P L {y)){P-P R {y)) 

The two-loop correlator is obtained immediately from (|A.13| ) 

oo 

W(L, L') = e 5 ( L+L ') nI n (2RL) I n (2RL') (A.21) 

n=l 

Its Laplace image depends only on the positions of the endpoints of the cut and is given 
by 

W(P ,P> )= VT^TO^i^g^ " 2 , (A.22) 

In the continuum limit 

L = £/a, P = P R _+az, P R {y)=P R -af(y), /(0) = M (A.23) 
eqs. ( |A.16|) and (|A.17|) imply 

u'(-*)=V (A-24) 



/(„) v/a(z + M) 

where u(z) = U(P R + az). If the potential u(z) scales as a a , the dependence on the cutoff 
a can be elliminated by rescaling u(z) — > a a u(z),y — > a a ~ l / 2 y. Eq. ( |A.24| ) is then the 



Laplace image of ( |4.14|) for the potential 

Finally, the two-loop correlator in the scaling limit reads 

-M(£+f) 

w{£, £') = VI e + £ , VF (A.25) 

= ^ = "7T 7T7 Mv^+M + v/TTM) (A.26) 

4 (z — z'Y az az' 

In the case of gaussian potential U($) = -$ 2 / 2 the one-loop amplitude can be cal- 
culated explicitely in all orders in k in the continuum limit |)5| |)6) 



-M 

w(£) = I dy(y\e l ^^ d / d y^\y) = r */\-MM\ (A27) 

-oo 
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